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1  Introduction 


The  primary  motivation  for  this  research  is  the  interest  by  the  Air  Force  and  many  other  organizations 
in  developing  and  deploying  large,  precise,  lightweight,  space-based  antennas  and  optical  telescopes.  Large 
diameter,  optical  quality  membrane  reflectors  may  well  be  the  critical  components  that  make  such  structures 
possible.  The  prevailing  paradigm  for  creating  doubly-curved  membrane  surfaces  is  the  pressurized  lenticular 
configuration,  illustrated  in  Figure  1.  However,  there  are  serious  difficulties  that  must  be  overcome  in  order 
for  such  a  configuration  to  be  successfully  deployed  in  space.  These  problems  are  discussed  at  some  length 
in  Reference  [1], 


Figure  1:  Optical  lenticular  for  imaging. 


An  alternative  to  the  pressurized  lenticular  configuration  is  a  thin  laminate  shell  consisting  of  a  membrane 
and  a  dielectric  coating  commonly  used  in  optics  to  meet  a  high  reflectivity  requirement.  This  laminate  shell 
is  referred  to  as  a  stress-coated  net-shape  membrane  reflector.  The  net-shape  membrane  portion  of  the 
laminate  is  realized  by  manufacturing  a  polymer  to  nearly  its  final  shape,  which  is  typically  either  spherical 
or  parabolic.  The  net-shape  process  involves  several  basic  steps.  The  solvent-based  polymer  is  initially  cast 
on  a  mandrel  (also  referred  to  as  a  “mold”)  having  the  desired  surface  qualities  (e.g.,  optically  smooth,  and 
having  the  required  shape),  where  the  solvent  is  allowed  to  evaporate,  leaving  a  thin  high  quality  membrane 
of  the  desired  shape.  At  this  point,  however,  significant  internal  stress  has  developed  in  the  membrane  due 
to  the  solvent  evaporation  process.  The  membrane  is  then  taken  through  an  annealing  process  in  which  the 
mold  and  membrane  are  heated  to  nearly  the  glass  transition  temperature  of  the  polymer,  eliminating  most 
of  the  shrinkage  stress.  However,  the  coefficient  of  thermal  expansion  (CTE)  of  the  membrane  is,  in  the 
systems  we  consider,  higher  than  that  of  the  mold.  As  the  system  is  allowed  to  cool  to  room  temperature, 
the  membrane  attempts  to  contract  more  than  the  mold,  due  to  its  higher  CTE.  Since  the  membrane  is  fully 
constrained  by  the  mold,  hence  cannot  complete  its  contraction,  the  effect  of  the  CTE  mismatch  is  to  induce 
a  (tensile)  thermal  stress  in  the  membrane.  Such  a  stress  is  referred  to  as  non-mechanical  (or  inelastic  or 
residual ),  i.e.,  it  is  a  stress  that  exists  in  the  absence  of  any  displacement-related  strain  (see,  for  example, 
Fung  [2],  pp.  354-355).  This  residual  thermal  stress  would  act  to  deform  the  membrane  from  its  initial  shape 
upon  removal  from  the  mold.  The  other  serious  problem  is  the  simple  fact  that  a  membrane  is  ” flimsy”,  that 
is,  it  lacks  the  stiffness  required  to  resist  bending  due  to  external  loads.  Examples  of  such  loads  are  gravity 
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and  wind  in  a  near  earth  environment,  and  slewing  or  other  forces  used  for  control  in  a  space  environment. 
At  any  rate,  upon  removal  from  the  mold,  a  net-shape  membrane  would  not  be  expected  to  retain  the  shape 
of  the  mold. 

The  research  reported  here  addresses  the  possibility  of  solving  both  the  CTE  mismatch  problem,  and 
lack  of  stiffness,  by  applying  to  the  membrane  a  coating  with  an  intrinsic  compressive  stress  designed  to 
compensate  the  CTE  mismatch  stress,  as  well  as  provide  enough  stiffness  to  maintain  the  desired  shape 
under  various  loads.  Specifically,  we  examine  the  effects  of  gravity  and  uniform  pressure  loads  on  the 
stress-coated  membrane.  In  Volume  I  of  this  two-volume  report  the  method  of  asymptotic  expansions  is 
used  to  derive  various  theories  of  stress-coated  membranes  from  the  general,  geometrically  nonlinear,  three- 
dimensional  theory  of  elasticity.  In  Volume  II  we  present  solutions  of  the  equations  satisfying  various  types 
of  boundary  conditions. 

2  Reference  Placement  and  Reference  Configuration 

We  introduce  a  region  C  of  3-dimensional  Euclidean  space  in  the  form  of  a  thin  right  circular  cylinder  of 
radius  a  and  uniform  thickness  (or  height)  h  «  o,  and  refer  to  this  purely  mathematical  construct  as  the 
reference  placement  of  a  coated  membrane  shell.  This  cylinder  is  further  assumed  to  be  divided  into  two 
coaxial  cylinders  of  the  same  radius  a,  one  of  thickness  /iff,  the  other  of  thickness  hc ,  so  that  h  =  h8  +  hc, 
as  shown  in  the  lower  portions  of  Figures  2  and  3.  We  assume  given  a  fixed  orthonormal  Cartesian  basis 
{i,j,k}  with  origin  O  at  the  center  of  the  circular  disk  defined  by  the  intersection  of  C  and  a  bisecting 
plane  orthogonal  to  the  axis.  An  arbitrary  point  P  of  the  reference  placement  may  be  specified  by  either  its 
Cartesian  coordinates  XA  =  {X,y,Z},  or  its  cylindrical  coordinates  QA  =  {R,  0,Z}.  Thus,  the  bisecting 
or  middle  plane  of  C  is  defined  by  Z  =  0,  and  the  axis  of  C  by  the  line  X  =  Y  =  0  through  O.  The  position 
vector  of  P  with  respect  to  O  is  given  by 

X  =  Xi  +  Yj  +  Zk  =  R  cos0i  +  R  sin©  j  +  Z k.  (2.1) 

We  introduce  orthonormal  basis  vectors  {Eh,E©,E£}  associated  with  the  cylindrical  coordinates,  defined 

by 

Ea=X,a/\Xa\,  where  X,  A  =  ^,  (2.2) 

so  that 

Eh  =  cos©i  +  sin©j,  E©  =  —  sin0i  +  cos0j,  E  z  =  k,  (2.3) 

in  terms  of  which  we  can  write  the  position  vector  as 

X  =  REr  +  ZEz.  (2.4) 

The  physical  system  of  interest  is  a  laminate  material  body  in  the  form  of  an  initially  curved  membrane 
substrate  to  which  an  optical  coating  has  been  applied.  The  reference  configuration  of  this  coated  membrane 
is  assumed  to  be  a  region  S  defined  by  a  mapping  <j)  from  the  reference  placement  C,  under  which  a  point  P 
of  C  is  mapped  to  some  material  point  P  =  <£(P)  of  5.  A  point  S  of  the  middle  plane  of  C  with  coordinates 
(it!,  0, 0)  is  mapped  by  <f>  to  a  point  S  of  the  middle  surface  of  S  with  coordinates  (i?,  0,  Zs)>  where 

ZS  =  r  (R),  (2.5) 

hence  we  are  assuming  that  the  middle  surface  is  a  surface  of  revolution.  The  azimuthal  coordinate  0  of 
any  point  of  C  is  assumed  unchanged  by  this  mapping,  so  that  0  =  0  on  5.  This  action  of  (j>  on  the  middle 
plane  is  illustrated  in  the  upper  portions  of  Figures  2  and  3. 

The  action  of  <f>  on  points  off  the  middle  plane  depends  on  the  distribution  of  the  thicknesses  hc  and  h8  as 
a  result  of  the  processes  used  to  cast  the  membrane  on  the  mold,  and  to  apply  the  coating  to  the  membrane. 
The  simplest  model  results  by  assuming  the  coated  membrane  to  have  constant  axial  thicknesses,  so  that 
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Surface  of  Mold 


Figure  2:  Definition  of  the  reference  configuration  S  (upper  part  of  Figure)  of  a  coated  membrane  shell  of  revolution 
as  a  mapping  from  the  reference  placement  C  (lower  part  of  Figure),  assuming  the  thicknesses  hc  and  hs  to  be  constant 
along  any  line  parallel  to  the  axis. 


an  arbitrary  point  P  with  non-zero  axial  coordinate  Z  is  mapped  to  a  material  point  P  with  coordinates 
{i?,  0,  Zj,  obtained  by  translating  a  distance  Z  (equal  to  the  original  axial  coordinate)  from  S  along  the 
axial  direction  E z,  as  shown  in  the  upper  portion  of  Figure  2.  The  complete  mapping  can  be  determined 
geometrically  as  follows.  The  position  vector  of  P  with  respect  to  O  can,  according  to  Figure  2,  be  written 
in  two  different  ways: 

X  =  REr  +  ZEz  =  Y  +  ZEz ,  (2.6) 

where  Y  is  the  position  vector  of  S  with  respect  to  O.  Now,  Y  ~  REr  +  T(R)  Ez,  hence  we  have  from 


3 


(2.6): 


REr  +  ZEz  =  REr  +  \T(R)  +  Z]Ez, 

and  a  comparison  of  components  on  both  sides  of  the  second  equality  yields  the  remaining  two  component 
mappings  of  <p  (0  =  ©  is  the  other  component  mapping): 

R  =  R,  and  Z  =  T(R)  +  Z.  (2.7) 

The  position  vector  of  a  material  point  in  the  reference  configuration  S  of  the  coated  membrane  is  thus  given 
in  terms  of  coordinates  on  the  reference  placement  C  by 

X  =  REr  +  {T{R)  +  Z]Ez.  (2-8) 

For  comparisons  with  finite  element  analyses  that  use  shell  elements  in  their  formulation,  it  is  perhaps 
more  appropriate  to  assume  the  coated  membrane  to  have  constant  thicknesses  hc  and  he  normal  to  the 
middle  surface.  In  this  case^an  arbitrary  point  P  with  non-zero  axial  coordinate  Z  is  mapped  to  a  material 
point.  P  with  coordinates  {R,@,Z},  obtained  by  translating  the  distance  Z  along  the  unit  normal  EN  to  the 
middle  surface  at  S  (see  Figure  3).  The  position  vector  of  P  with  respect  to  O  is  then  given  by 

X  =  REr  +  ZEz  =  Y  +  ZEN,  (2.9) 

where  Y  =  REr  +  T(R)  Ez,  as  before.  To  compute  the  unit  normal  Ejv,  we  first  note  that  equation  (2.5) 
can  be  written  as  ip{R,  Zs)  =  0,  where  ip  is  the  function  defined  by 

iP(R,  Z)  =  Z  -  T(R).  (2.10) 


The  unit  normal  to  the  midsurface  is  the  normalized  gradient  of  this  function: 

Vip  —T'REr  +  Ez  mn 

E"  -  W\  ~  v/i  +  <r,*)’  ’  (  ’ 

where  here  is  the  ordinary  derivative  of  T  with  respect  to  R  (note  that  this  slope  is  negative  in  the  first 
quadrant).  Substituting  these  results  in  (2.9)  yields 


REr  +  ZEz  =  REr  +  T(R)EZ  +  Z 


W1  +  (r  *F  )  5 


from  which  the  component  mappings  of  <j>  are  given  by 

Ja _ 


i?  =  R-  Z 


=  R  +  Z  sin  a, 


Vi  +  M* 

z  =  T(R)  +  Z  y=  =L==  =  T(R)  +  Z  coso, 


(2.12) 


Vi  +  (r,fl)2 

where  a  is  the  angle  between  Ez  and  the  unit  normal  Ejv  at  a  point  on  the  middle  surface  (cos a  =  Ez-Ejv). 
From  (2.12)  we  have 


2  P  '  j  . . 

cos  a  =  — ,  sin  a  =  — where  $  =  y  1  +  (r,j{)2  . 


(2.13) 


The  position  vector  of  a  material  point  in  the  reference  configuration  of  the  coated  membrane  is  thus  given 
in  terms  of  coordinates  on  the  reference  placement  by 

X  =  -  .Z^jp)  Efl  +  (r(R)  +  Z^j  Ez  =  (R  +  Z  sin  a)  E^  +  [T(R)  +  Zcosa]Ez.  (2.14) 


4 


Surface  of  Mold 


Figure  3,  Definition  of  the  reference  configuration  S  as  a  mapping  from  the  reference  placement  C,  assuming  the 
thicknesses  hc  and  hs  to  be  constant  along  any  line  through  a  normal  to  the  middle  surface  Zs  =  r(J?). 


The  use  of  a  model  with  constant  thicknesses  normal  to  the  middle  surface  is  considerably  more  compli¬ 
cated  to  analyze  theoretically  than  one  with  constant  axial  thicknesses.  In  the  remainder  of  this  work  we 
use  the  simpler  model.  To  get  an  estimate  of  the  errors  made  in  choosing  the  constant-axial-thickness  model 
over  the  constant-normal-thickness  model,  we  consider  the  important  case  of  a  paraboloidal  middle  surface 
defined  by 


T(R)  =  To  -  ±R2, 


(2.15) 


where  /  is  the  focal  length  and  T0  is  the  apex  displacement.  Since  T(a)  =  0,  it  follows  that  T0  =  a2/(4/), 
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hence 


(2.16) 


m-  j]  (<?-#)■ 

The  tangent  of  the  slope  angle  is,  from  (2.13),  tana  =  — r  r,  so  that 

tan  a  =  R. 

The  /-number  of  the  paraboloid,  which  we  denote  by  F&,  is  defined  by 

F*  =  — 

2a’ 

so  we  can  write  the  last  two  equations  in  terms  of  /-number  as 

T{R)  =  8^F*  (°2  "  ^  ’  tanQ  =  ~T’R  =  4 ^F*K 
The  angle  a  has  its  maximum  value  at  the  edge  R  =  a,  in  which  case 

1 

tan  Olmax  —  ' 

For  the  optical  applications  envisioned  here  we  expect  to  have  /-numbers  of  2  or  greater,  hence 


tan< 


(2.17) 

(2.18) 

(2.19) 

(2.20) 

(2.21) 


Use  of  the  constant-axial-thickness  model  entails  approximating  cos  a  «  1  and  sin  a  ps  0,  so  for  this  lowest  /- 
number  one  might,  expect  any  differences  between  our  theoretical  results,  and  finite  element  analysis  results, 
to  be  at  least  partly  attributable  to  these  approximations. 


3  Deformation,  Displacement,  and  Strain 

When  the  coated  membrane  shell  is  removed  from  the  mold,  releasing  it  from  its  constraints,  it  deforms  until 
a  new  equilibrium  configuration,  which  we  refer  to  as  the  deformed,  or  current,  configuration,  is  attained. 
The  deformation  is  assumedto  be  described  mathematically  by  a  one-to-one  invertible  mapping  /  that  maps 
the  body  point  located  at  P  of  the  reference  configuration  C  to  a  new  point  p  =  /( P).  The  set  of  image 
points  of  /  defines  the  deformed  configuration  T>.  Variables  that  refer  to  points  of  this  new  configuration 
will  be  denoted  by  lower  case  Latin  letters,  e.g.,  xa  =  {x1  ,  x2 ,  x3  }  =  {x,  y,  z}  are  Cartesian  coordinates,  and 
qa  =  {g1,  g2,  q'3}  =  {r,  0,  z)  are  cylindrical  coordinates  on  the  deformed  configuration.  Thus,  in  terms  of 
these  cylindrical  coordinates,  the  mapping  p  =  /( P)  is  coordinatized  by 

9°(p)  =  q°(f( P))  =  (9a°/)( P)  s  /“( P),  (3.1) 

where  the  functional  compositions  qa  o  /  =  /“  define  the  cylindrical  component  mappings  of  the^ mapping  /. 
Assuming  the  arbitrary  point  P  to  be  coordinatized  by  the  cylindrical  coordinates  QA  =  {R, ©,  Z}  on  the 
reference  configuration,  equation  (3.1)  can  be  written  as 

9°(p)  =  fa[QHP),Q2(HQHn-  (3-2) 

It  is  more  convenient,  however,  to  relate  points  of  the  deformed  configuration  to^  points  of  the  reference 
placement  S.  Each  point  P  of  C  is  the  image  of  some  point  P  of  S  via  the  mapping  P  =  </>( P).  This  mapping 
can  be  used  to  express  the  actual  deformation  of  the  shell  in  terms  of  points  of  the  reference  placement,  viz., 

P  =  /( P)  =  7(*(P»  =  (7°  M  P)  =  /( P),  (3-3) 


6 


where  the  functional  composition  /  =  /  o  (j>  defines  a  mapping  from  the  reference  placement  to  the  current 
configuration.  The  three  mappings  we  have  introduced  are  illustrated  in  Figure  4.  Thus,  cylindrical  coordi¬ 
nates  of  points  on  the  current  configuration  are  given  in  terms  of  the  cylindrical  coordinates  QA  —  {i?,  0,  Z} 
on  the  reference  placement  by 

3°(p)  =  /a[Q1(P),Q2(P),Q3(P)],  (3.4) 


analogous  to  (3.2). 


Figure  4:  Mappings  relating  the  reference  placement,  reference  configuration,  and  deformed  configuration. 

The  position  vector  of  a  material  point  in  the  reference  configuration  S  of  the  coated  membrane  is  given 


in  terms  of  coordinates  on  the  reference  placement  C  by  equation  (2.8),  repeated  here: 

X  =  REr  +  [F(i?)  +  Z]  Ez.  (3.5) 

From  (2.3),  also  repeated  here: 

Eh  =  cos0i  4-  sin©j,  E©  =  —  sin0i  4-  cos0j,  E z  =  k,  (3.6) 

we  obtain  the  differentials  of  the  basis  vectors,  expressed  in  terms  of  the  same  basis  vectors: 

dEtR  =  d0E©  5  dE©  =  —  d©E_R,  dE^  =  0.  (3.7) 

Using  these,  we  find  for  the  differential  of  the  position  vector: 

dX  =  diiEjj  4-  i?d©E@  4-  [dZ  4-  T^dii]  Ez*  (3.8) 

In  (3.8)  we  introduce  the  differential  forms  Qr  =  dR ,  0©  —  RdQ,  and  £lz  =  dZ ,  to  write  it  as 

dX  =  Oh  Eh  4-  0©E©  4*  [Qz  ^~E^rQr]Ez  =  O^E^  4-  F^QhE^,  (3.9) 


where  the  usual  summation  convention  on  repeated  indices  is  to  be  understood,  as  in  the  first  term  of 
the  second  equality  of  (3.9),  unless  otherwise  stated.  Thus,  the  differential  forms  0^  can  be  written  as 
Qa  =  HAdQA  (no  sum  on  A),  where  Hr  =  1,  H©  =  R ,  and  Hz  —  1  are  referred  to  as  scale  factors. 
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The  position  vector  of  a  point  on  the  deformed  configuration  is  given  in  terms  of  its  cylindrical  coordinates 

by 

x  =  r  cos0i  +  r  sin0j  +  z k  =  rer  +  zez,  (3.10) 

where 


er  =  cos0i  +  sin0j,  e$  =  —  sin0i  +  cos0j,  ez  =  k.  (3.11) 

Similarly  to  (3.9),  the  differential  of  this  position  vector  can  be  written  as 

dx  =  uaea,  (3.12) 

where  uia  =  hadqa  (no  sum  on  o)  sure  differential  forms  on  the  deformed  configuration  (with  hT  =  hz  =  1, 
he  =  r). 

Prom  equation  (3.4),  we  obtain 


dqa 


i 


or,  replacing  coordinate  differentials  by  their  respective  differential  forms, 


hg  0fa 
HAdQA 


(3.13) 


where 


FaA  = 


hg  df* 
Ha  9Qa 


(no  sum  on  either  a  or  A), 


(3.14) 


are  the  elements  of  the  matrix  F  of  deformation  gradients  from  the  differential  forms  on  the  reference 
placement  to  differential  forms  ^n  the  deformed  configuration. 

The  displacement  of  point  P  of  the  reference  configuration  to  point  p  of  the  deformed  configuration  is 
defined  by  the  vector  field 


U  EE  X  -  X. 


Prom  this  relation  we  obtain 


(3.15) 


dx  =  dX  4  du.  (3.16) 

If  we  write  the  displacement  field  u  in  terms  of  its  components  Ua  in  the  orthonormal  cylindrical  basis  of 
the  reference  placement,  i.e., 


u  =  UrEr  4  t/©  E©  4  UzEz ,  (3.17) 

where  the  components  Ua  are  assumed  to  be  functions  of  the  reference  placement  coordinates  {ii,  0,  Z},  we 
find  for  its  differential: 


du  =  (UR,RdR  4  J7R,©d0  4  Ur,z&Z )  "Er  4  Ur  (IEr 
4  ( Ue.RdR  4  J7©,©d0  4*  Us.zdZ)  E©  4  Ue  dE© 
4  (Uz,RdR  4  Uz,ed&  4  Uz,zdZ)  Ez, 

=  UR,RQR  +  (UR,ex  — )  +  Ur,zCIz  er 


4 

4 


Ue,R^R  4 


^e,e  +  tfo) 


ft©  4  Ue.z^z 


E© 


Uz’rQr + (^) 


fl©  +  Uz,z&z 


Ez, 


(3.18) 
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where  U a, b  denotes  the  partial  derivative  of  Ua  with  respect  to  Qb-  This  can  be  written,  using  the 
summation  convention,  as 


du  =  HabQbEa, 


(3.19) 


where  Hab  —  Ua-,b  are  the  elements  of  the  matrix  H  containing  the  components  Ua,b  of  the  covariant 
derivative  of  u,  distinguished  by  a  semicolon  in  place  of  the  comma.  Explicitly,  we  have 


Ur,r  {Ur,q-Uq)/R  Ur,z 
H  —  VetR  (Ue,e  +  Ur)/R  Uq >z 
Uz,r  Uz,e/R  Uz,z 


(3.20) 


The  Green-Lagrange  strain  tensor  E  is  defined  by 


®  —  2  '  dx  dX  •  dX^  —  —  ^dX  ■  du  +  du  •  dX  +  du  •  du^  ,  (3.21) 

where  (3.16)  was  used  to  get  the  second  equality.  From  (3.9)  and  (3.19),  we  have 

dX-du  =  (ficEc  +  T,rSIr-Ez).HabSIbVa  =  {HAb +T,rHzb5ar)^a^b, 
after  relabeling  of  dummy  summation  indices,  and  similarly, 

du-dX  =  Hab^bEa  •  (Ac  Ec  +  rHfiflE2)  =  (Hba  +  T'RHzaSbrWaQb. 

Substitution  of  the  last  two  results  in  (3.21)  yields 

®  =  2  -  ^AB  +  ^ba  +  T'RHzbSar  +  TtRHZA$BR  +  HcaHcb  )  =  EabQa&b,  (3.22) 

where  Eab  =  ( Hab  +  Hr  a  +  T ,rHZb$ar  +  T  ,rHZa$br  +  HcaHcb) /2  are  the  elements  of  the  Green- 

Lagrange  strain  matrix  E.  Carrying  out  the  algebra,  we  obtain  the  following  expressions  for  the  components 
of  the  strain  tensor  in  cylindrical  coordinates  on  the  reference  placement: 

Err  =  UR.R  +  T,rUz,r  +  ^  (f/^%  +  U^R  +  U^R  )  ,  (3.23) 


_  k©,e  +  Ur  1  ( ©  ~  kr©  )2  +  ( Ue,e  +  Ur  )2  +  U%  B 


Ezz  =  Uz,z  +  \  (Ur  Z  +  uiz  +  JJZ  Z  )  , 


Er&  -  -  Ub,r  +  r ~’9  +  UR,e  —  —  +  Ur’r(Ur&  ~  u&)  +  ^e,fl(lde,e  +  Ur)  +  Uz,'r  Uz^e 


Eez  =  1  U@  z  +  ^ R*z(U r,q  —  Ue)  +  Ue,z(Uete  +  Ur)  +  Uztz  Uz,s 


Erz  —  2(URiZ  +  Uz,R  +  T,rUz,z  +  Ur^rUr^z  +  U@,RUBiz  +  Uz,RUz,z)-> 


9 


where  Eqr  =  Erq ,  Eze  =  Eez,  and  Ezb  =  Erz,  i.e.,  E  is  symmetric.  It  is  important  to  note  that  we 
are  here  taking  over  from  the  classical  theory  of  laminates  the  fundamental  assumption  that  the  displace¬ 
ment  components,  hence  the  strain  tensor  components,  are  continuous  through  the  coated  membrane  shell 
laminate. 

Using  (3.12),  (3.13),  (3.9),  and  (3.19)  in  (3.16),  we  obtain 

FaB^B^a  =  ft  A  Ea  +  T^rQrEz  4“  HaB^B^a  =  {$AB  +  T^rSaZ^BR  +  Hab)&bEa, 
from  which  follows  the  useful  relation: 

e0  =  (<Sab  4-  T^rSazSbr  4*  Hab  )  (F~1)bo.Ea  =  Kab{F~1)boFia  =  C^Ea,  (3.29) 

where 

Kab  =  Sab  4*  r  rSazSbr  4-  Hab ,  Oa<x  =  Ea  -ea  =  Kab(F  1)bc-  (3.30) 

The  matrix  O  with  elements  defined  by  (3.30)  must  be  orthogonal,  satisfying  0T0  =  00T  =  I  (where 
the  T-superscript  denotes  a  transposed  matrix),  since  both  bases  {ea}  and  {Ea}  are  orthonormal.  It  is  an 
example  of  a  shifter  [3,  p.  9],  in  this  case  from  one  orthonormal  basis  to  another.  Assuming  both  bases  to 
have  been  chosen  as  right-handed,  the  determinant  of  O  must  be  1,  i.e.,  det(O)  =  det(0T)  =  1.  From  the 
matrix  form  of  the  second  equation  of  (3.30),  i.e.,  O  =  KF~l  =>  F  =  0T if,  it  then  follows  that 

J  =  det(F)  =  det(Jf),  (3.31) 


where  J  =  det (F)  is  the  Jacobian  determinant  of  the  matrix  of  deformation  gradients.  It  is  easy  to  show 
from  (3.6)  and  (3.11)  that  the  shifter  from  the  orthonormal  cylindrical  basis  on  the  reference  placement  to 
the  orthonormal  cylindrical  basis  on  the  current  configuration  is  given  by 


O  = 


cos  (0  —  0)  sin  (9  —  0) 
—  sin  (6  -  0)  cos  (0  —  0) 
0  0 


0" 

0 

1 


(3.32) 


Since  F  =  0TK  we  obtain,  using  (3.30)  and  (3.32): 


'cos  (0  -  0) 

-  sin  (0  —  0) 

O' 

1  +  Ur,r 

(UR,e  -  Ue)/R 

UR,z 

F  = 

sin  (0  -  0) 

cos  (0  -  0) 

0 

Ue,R 

1  +  (Ue,e  +  Ur)/R 

Ue,z 

0 

0 

1 

T,r  +  Uz,r 

Uz,e/R 

1  +  Uz,zm 

which  yields 

FrR  =  (1  4*  Ur,r)  cos  (0  -  0)  -  Ue,R  sin  (0  -  0),  (3.34) 


Fre  =  (Urt'eR  — )  cos  (0  -  0)  -  (l  +  t/e,eR+— )  sin  («-©).  (3.35) 

Frz  =  UR,z  cos  (0  -  0)  -  Ue,z  sin  (0-0),  (3.36) 

For  =  (1  +  Ur,r)  sin  (0  -  0)  +  Uq,r  cos  (0-0),  (3.37) 

Foe  =  (g”g_zge)  sin  (0  -  0)  +  (l  +  Ue'e  +  ^)  cos  (0  -  0),  (3.38) 


Fez  =  Ur,z  sin  (6  —  0)  4-  Us,z  cos  (6  —  0), 


(3.39) 


U7  A 

=  f>  +  Fze  —  ,  Fzz  —  1  +  Uz,z*  (3.40) 

Note  that  dv  =  JdF,  where  dV  is  a  volume  element  on  the  reference  placement  and  dv  is  a  volume  element 
on  the  current  configuration.  In  order  for  the  matrix  F  of  deformation  gradients  to  be  invertible,  we  disallow 
the  possibility  that  J  =  0.  We  also  disallow  the  possibility  that  J  <  0,  as  that  would  imply  that  a  volume 
element  could  have  a  negative  volume.  Thus,  we  require  that 


where,  from  (3.30)  and  (3.31): 


j  =  i  +  Ur,r  +  uz,z  +  u°-e  +  u*  +  (gas-ps)  Vzz .  ^f-ue,z  +  UreUZ:Z 

+  (U°*  +  U*')  -  ^f-Ue,zUR,R 

-  -  (^)  ♦  (2^)  r, 

+  (— £fe,z  UXtR  -  Un,z r,R  -  tfz,*l/*,z  -  (t/e'eR+^  [/r,z  rfl 

-  ([7e  ejf^il)  +  ^UR>zUe,R. 


4  Equilibrium  Equations 

Equilibrium  of  a  deformed  body  requires  that  both  the  net  force  and  net  moment  of  force  on  any  part  of  the 
body  vanish.  The  vanishing  of  the  net  moment  is  well-known  to  imply  the  symmetry  of  the  Cauchy  stress 
tensor  er  (our  notation  for  stress  tensors  follows  that  of  References  [4,  pp.  134-136],  and  [5,  Chapter  4]). 
Here,  we  begin  with  the  force  equilibrium  equations  written  in  terms  of  the  Cauchy  stress  tensor,  and  then 
reformulate  them  in  terms  of  the  first  and  second  Piola-Kirchhoff  stress  tensors.  Let  V  denote  the  volume 
of  an  arbitrary  part  of  the  coated  membrane  in  its  plate-like  reference  placement ,  and  denote  by  dV  the 
boundary  surface  of  this  part.  Under  the  deformation  /  defined  by  (3.3),  V  is  mapped  to  f(V ),  bounded 
by  the  surface  df{V).  In  the  presence  of  a  gravitational  body  force  f5,  force  equilibrium  of  the  arbitrary 
deformed  volume  f{V)  requires  that 

f  cr  *n da  +  /  f gdv  =  <b  oabnteada  A  /  pgdv  =  0,  (4.1) 

Jdf(V)  Jf{v)  Jdf(v)  Jf(v) 

where  er  =  oabeaebj  n  =  ncec  is  the  unit  normal  to  the  deformed  surface,  p  is  the  mass  density  of  the 
material,  and  g  =  gez  =  gEz  is  the  gravitational  acceleration  (assumed  to  act  in  the  “up55  direction  along 
the  positive  Z-axis  in  Figure  2),  expressed  in  the  orthonormal  cylindrical  bases.  The  surface  integral  can 
be  written  in  terms  of  quantities  on  the  reference  placement  using  a  version  of  Nanson’s  formula  (see,  for 
example,  [6],  p.  249,  or  [7],  p.  88),  viz., 

nada  =  JiF-^AaNAdA,  (4.2) 

where  the  Na  are  components  in  the  orthonormai  cylindrical  basis  of  the  unit  normal  N  to  the  surface 
element  of  area  dA  in  the  reference  placement.  The  volume  integral  is  transformed  to  one  over  the  reference 
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placement  volume  by  the  substitution  dv  -  JdV,  where  dV  is  the  reference  placement  volume  element. 
Making  these  substitutions  in  (4.1)  yields 


<p  CT0j,ni,eada+  f  pgezdv  =  ®  a  abJ(F  1)AbNAeadA  +  /  JpgezdV 
Jdf(V)  JfCP)  J&V  -If  (?) 

=  <fc  PaANAeadA  +  f  Pagezdv  -  0,  (4.3) 
JdV  JfCP) 


where  po  =  Jp  is  the  mass  density  of  material  in  the  reference  configuration,  and 

PaA  =  J<7ab(F~1)Ab, 


(4.4) 


are  the  components  of  the  nonsymmetric  first  Piola-Kirchhoff  stress  tensor  P,  which  appears  naturally  in 
transforming  from  deformed  configuration  surface  elements  to  reference  placement  surface  elements.  How¬ 
ever,  (4.3)  is  expressed  in  term  of  components  along  the  deformed  configuration  basis  vectors  eQ.  The  second 
Piola-Kirchhoff  stress  tensor  S  arises  naturally  by  using  (3.29)  to  shift  to  components  along  the  orthonormal 
cylindrical  basis  vectors  Eb  of  the  reference  placement,  obtaining  for  the  surface  integral  in  (4.3): 


<f  PaANAeadA  =  <f  PaANAKBc{F  1)caEBdA  =  ®  KBcScA^A^BdA,  (4.5) 

JdV  JdV  J  9V 


where  we  identify 

ScA  =  {F^CaPaA  =  J  (F-')Ca<Tab{F-l)Ab  (4-6) 

as  the  components  of  the  second  Piola-Kirchhoff  stress  tensor.  Thus,  we  can  write  (4.5)  as 

<£  PaANAeadA  =  <f  KBcScANAEBdA  =  <£  TBANAEsdA,  (4.7) 

JdV  JdV  Jdv 

where  it  was  convenient  to  introduce  yet  another  (nonsymmetric)  stress  tensor  T  with  components  defined 
by 


Tba  =  KbcSca  =  (£bc  +  r  ,r$bz8cr  +  Hbc)Sca- 


(4.8) 


Using  (4.7)  in  the  force  equilibrium  equations  (4.3)  yields 

{  TBANAEBdA  +  f  p0gEzdV  =  0.  (4.9) 

Jdv  JfCP) 

From  (4.9)  we  infer  that,  just  as  the  equations  of  equilibrium  in  terms  of  the  Cauchy  stress  follow  by  an 
application  of  the  divergence  theorem  to  convert  the  surface  integral  on  the  deformed  configuration  to  a 
volume  integral,  yielding  from  (4.1)  local  equilibrium  equations  of  the  form 

(Tab;b  +  pg^az  =  0,  (4-10) 

where  the  left-hand  side  includes  the  covariant  divergence  of  a  in  cylindrical  coordinates  on  the  deformed 
configuration,  so  also  must  the  nonsymmetric  stress  tensor  T  satisfy  local  equilibrium  equations  of  the  form 

Tab-,b  +  Po9Saz  —  0,  (4.11) 

where  the  left-hand  side  includes  the  covariant  divergence  of  T  in  cylindrical  coordinates  on  the  reference 
placement.  Taking  into  account  the  nonsymmetric  nature  of  the  components  TAb,  the  component  equations 
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5  Constitutive  Relations 


Prior  to  applying  a  coating,  a  non-mechanical  thermal  stress  due  to  CTE  mismatch  between  mold  and 
membrane  develops  as  the  two  cool  to  room  temperature  following  the  annealing  process.  We  denote  this 
mismatch  stress  by  A  coating  of  thickness  hc  is  then  applied  to  the  constrained  membrane,  and  assumed 
to  be  perfectly  bonded  to  it,  hence  fully  constrained  geometrically.  The  coating  typically  undergoes  some 
microstructural  change  during  the  coating  process,  inducing  an  intrinsic  coating  stress  S™ m  (which  may  be 
either  tensile  or  compressive).  If  the  membrane  is  coated  at  a  temperature  different  from  the  temperature 
at  which  the  mold  is  stress-free,  a  thermal  mismatch  stress  may  also  appear  in  the  coating. 

We  assume  that  both  materials  are  linearly  elastic,  uniform,  homogeneous,  and  isotropic,  and  that  they  re¬ 
main  perfectly  bonded  after  removal  from  the  mold  (guaranteeing  continuity  of  the  displacement  components 
across  their  interface).  As  in  [9],  [10]  and  [11],  we  have  begun  with  the  full  three-dimensional  Green-Lagrange 
strain  tensor  Eab>  including  all  geometrically  nonlinear  terms  involving  displacement  components  or  their 
partial  derivatives,  as  shown  in  equations  (3.23)-(3.28).  The  choice  of  the  Green-Lagrange  strain  tensor 
dictates  a  material  rather  than  spatial  description  of  the  deformation,  and  we  assume  that  in  each  of  the 
two  materials  we  have  a  simple  uniform,  linear,  isotropic  constitutive  relation  between  the  Green-Lagrange 
strain  tensor  and  the  second  Piola-Kirchhoff  stress  tensor: 

Sabi  =  S^Sab  +  [ Eab  +  (Eriu  +  Eeei  +  Ezzi)  $ab  j,  (5.1) 


where  the  subscript  i  denotes  the  ith  layer  of  the  laminate,  and  the  constants  Ei  and  Vi  are  Young’s  modulus 
and  Poisson’s  ratio,  respectively,  of  the  component  materials.  The  first  term  in  (5.1)  is  a  simple  way  of 
including  the  residual  or  non-mechanical  stress  in  material  i,  which  is  assumed  to  be  uniform  and  isotropic 
with  constant  value  Sfm.  It  is  a  slight  generalization  of  the  constitutive  relation  given  by  Fung  [2,  pp. 
354-355]  intended  to  account  for  a  thermal  stress,  in  which  case  it  is  given  by 


S?m  = 


Ej 

1  —  2  Vi 


nm 

> 


(5.2) 


where  efm  =  c*j  AT  is  the  thermal  strain,  a*  is  the  CTE  of  material  i,  and  AT  =  T  -  To  is  the  temperature 
deviation  from  some  reference  temperature  To-  It  is  assumed  here  to  include  possible  intrinsic  stresses  in 
the  materials  so  that,  in  general,  the  residual  stress  is  a  sum  of  intrinsic  and  thermal  stresses.  The 
component  forms  of  the  constitutive  relations  (5.1)  cam  be  written  as 


Snm  =  S?m  +  t 

-<  [  (1  -  Vi)  Err  +  Vi  ( Eee  +  Ezz )  ] , 

(5.3) 

Seei  =  S?m  +  t 

•«  [  (1  -  ^)  Eee  +  Vi  (  Err  +  Ezz )  ] , 

(5.4) 

Szzi  =  Sr  +  t 

\  [  (1  —  Vi)  Ezz  +  Vi  (  Err  +  Eee )  ]  i 

(5.5) 

Siiei  =  GiERe, 

Srz>  =  GiERZ,  Sezi  =  GiEez, 

(5.6) 

where  we  have  introduced  for  convenience: 


Si 


Ei 

(1 +  1/0(1-21/0’ 


Gi 


Ej 

1  +  ! >i 


(5.7) 


6  Boundary  Condition  of  Pressure 

To  formulate  appropriate  boundary  conditions  for  equations  (4.12)-(4.14),  we  consider  a  shell  that  is  in 
equilibrium  under  a  difference  in  hydrostatic  pressure  on  the  images  (under  the  deformation  /)  of  the  reference 
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placement  faces  Z  =  ft/2  and  Z  =  —ft/2.  Such  a  pressure  is  always  normal  to  the  deformed  surfaces,  and 
oppositely  directed  to  their  outward  unit  normal  vectors.  Thus,  traction  boundary  conditions  are  most  easily 
stated  in  terms  of  the  Cauchy  stress  tensor  on  the  deformed  surfaces: 

°'±*n±  =  -p±n±  =>  o±bnfda±e±  =  —  p±  nfda±e±,  (6,1) 

where  the  +  and  -  superscripts  denote  evaluations  either  on  the  planes  Z  =  h/2  and  Z  =  —ft/ 2,  respectively 
or,  as  in  (6,1),  on  their  images  under  /,  and  we  have  included  for  convenience  the  area  elements  da±  in  both 
sides  of  the  second  form  of  the  boundary  conditions.  Using  (3,29)  and  (4.2),  the  second  equation  of  (6,1) 
can  be  reformulated  in  terms  of  quantities  defined  on  the  reference  placement  as  follows: 

f0E±  =  -p±J±(F-1)^aNidA±K$B(F-^aE% 

which  reduces,  after  applying  the  definitions  (4.6)  and  (4.8),  to 

TcA*A  =  ~P±  (6.2) 

The  matrix  products  on  the  right-hand  side  of  (6.2),  viz.,  KF~1F~T,  where  F~T  is  the  transposed  inverse 
of  F,  are  easily  manipulated  to  the  identity  KF~~lF~T  =  K~T,  using  the  relation  O  =  KF~X  to  replace 
F~1 ,  and  the  orthogonality  of  the  shifter  O.  Thus,  after  a  relabeling  of  indices,  the  boundary  conditions 
take  the  form 

=  -p±J±(K-T)%Ni.  (6.3) 

However,  we  note  that  the  outward  unit  normals  to  the  faces  Z  —  h/2  and  Z  =  —ft/2  are  N+  =  E z  and 
N“  =  —Ezj  respectively,  so  that  =  ±Sbz  on  both  sides  of  (6.3),  hence  (6.3)  reduces  to 

T%z  =  -p±J±(K-T)±z.  (6.4) 

Recalling  that  J  =  det(F)  =  det(F),  it  follows  that  J±(iT_T):i:  is  just  the  matrix  of  cofactors  of  K±.  Thus, 
according  to  (6.4),  elements  of  the  third  columns  of  the  matrices  and  —p±J±(K~T)±  must  be  identical, 
yielding  the  final  forms  of  the  boundary  conditions  of  pressure: 


T»=t 

1 RZ 


T±  — 
1ez  — 


-P~ 


J7±  T7± 
uzieue,R 


R± 


-Uzjt 


1  + 


Ur  +  Uje 
R± 


-&r(1  +  uin)+ut* 

Ui  +  u$>e\ 


rp± 

Lzz 


(1+tfJb)  1+ 


JR± 


(6.5) 

(6.6) 

(6.7) 


7  On  the  Derivation  of  Theories  of  Two-Dimensional  Elastic  Bodies  from  the 
Three-Dimensional  Theory  of  Elasticity,  Using  the  Method  of  Asymptotic 
Expansions 

We  introduce  dimensionless  coordinates  (p,  £),  and  scaled  displacement  components  (U,V,W),  a  scaled 
reference  configuration  function  T,  scaled  stress  components  Sabi,  scaled  residual  in-plane  stresses  Si,  scaled 
pressure  loads  p ± ,  and  a  scaled  gravitational  body  force  poW-  defined  by  the  following  relations: 

R  =  ap,  Z  =  =  £oCi 

T  =  £rof ,  Uz  =  £maW,  Ur  =  elaU ,  U@  =  elaV, 
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(7.1) 

(7.2) 


.  SRRi  =  enE  iSRRU  Seei  =  e^Seeu  SRei  =  e^Sneu  (7.3) 

SRzi  =  epSt5^zi,  S©z*  =  epEtS©zi,  (7.4) 

Szz*  =  e9EjSzzi,  P±  =  ^E^j?*,  a/?0i  p  =  E<  e*  £,  (7.5) 

where  the  E<  are  two  arbitrary  constants  with  dimensions  of  stress  (E+  =  ESy  and  E“  =  Ec),  and  e  =  h/a  is 


the  (assumed  to  be  small)  scaling  parameter.  The  exponent  values  are  arbitrary  at  this  point,  but  typically 
satisfy  the  inequalities  r<m<t<n<t<p<q.  It  should  be  noted  that  the  original  variables  are 
functions  of  R ,  0,  and  Z,  e.g.,  UR  =  Ur(R>  0,  Z),  while  the  scaled  variables  are  all  functions  of  p,  0,  and  £, 
e.g.,  U  =  [7(p,  0,  C)  =  0,  hsC)*  From  these  definitions,  we  obtain  the  following  expressions  for 

the  partial  derivatives  of  the  displacement  components,  and  the  ordinary  derivative  of  the  surface-defining 
function  T: 

Ur,r  =  elU,p,  UR,e  =  6*aU,e,  UR,Z  =  (7.6) 

=  e%,  Ue,e  =  e‘aV,e,  Ue,z  =  (7.7) 

EA*,*  =  emW<p,  Uz,e  =  emaW,e,  Uz,z  =  em~lW,0  r,*  =  erftP.  (7.8) 

Substitution  of  these  expressions  into  equations  (3.23)-(3.28)  for  the  strain  components  yields: 


Err  =  e(U,p  +  £r+mf,pW,p  +  \  [e2"W2  +  £21  (d2p  +  V2)]  , 


Eee 


(V,e  +  U) 

w2 

+  e2m^4  +  e2* 

\u,e  -  V)2  +  (V,&  +  U)2' 

V  P  ) 

1  2  p> 

2p2 

Erq 


EZZ  =  £m~lW,<  +  e2m_2  \w$  +  e2t~2 J  (02<  +  V%), 
U,e-V\  +  £r+m  f,piy,e  ^ 


•sH'-*2*?2) 


+  e 


2/ 


+  £ 

P 

UAU,s-V)  +  VAV.e  +  U)' 

P 


(7.9) 

(7.10) 

(7.11) 


(7.12) 


Eez 


=  5 


(+£»Es+£— + 

P  P 


(7,c((7,e-Vr)  +  V.c(Ke  +  (7) 
P 


]}■ 


(7.13) 


Erz  =  ±  [e*-1  UA  +  e"W,P  +  +  e2™-1  +  e2'"1  ( )] .  (7.14) 

Substitution  of  (7.1)-(7.5)  in  (3.34)-(3.40)  yields  for  the  scaled  deformation  gradient  matrix  elements: 

FrR  =  (1  +  elU,p)  cos  (i 9-0 )  -  £%  sin  (9  -  0),  (7.15) 


Fr©  =  e*  — — —  )  cos  (0  -  0)  - 

V  p 


') 


sin  (0-0), 


(7.16) 


Frz  =  e*  1C7,^  cos  (0  -  0)  -  el  sin  ( 9  —  0), 


t-  if 


=  (1  +  elU,p)  sin  (9  -  0)  +  £lViP  cos  (9  -  0), 


(7.17) 

(7.18) 
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Fee  =  (  — — —  sin  (6  —  0)4-  1  +  el  (  ^?e  +  ^  cos  (0-0), 


Fgz  =  e*  1I/^  sin  (0  -  ©)  +  el  1 cos  (0-0), 


F,fl  =  er  T,fi  +  W,p>  Fz&  =  sm FzZ  =  1  +  e WiC. 


The  scaled  Jacobian  determinant  (3.42)  takes  the  form 


J  =  1  +  e”*-1  W*  +  e*  (t/p  +  V’e^+Lf  j  +  wx  \u,p  +  v'e^uj  -  wJj,  c  -  ^v;c 

+ ^  k  ( £sl±£)  -  v.  ( a^)l  e— ‘  w. 


+  £2,+”-‘  [  j  -  S'A)  - 

+  ^  {v„b,<  -  t?A)  I  +  =2,+r-1  f  „  f  ( 


17.  e  -  V 


(vPwx  -  wJt) 


+  (vA<  -  U,PV, c)  J  +  e2€+r-1  f  >p  j  _  ^^®±£  j  fjA  . 

The  scaled  expressions  for  the  strain  tensor  components  are  to  be  substituted  in  the  constitutive  relations 
(5.3)-(5.6).  We  begin  with  (5.5)  for  Szz%  from  which  we  obtain 


e*HiSzzi  =  enHiSfm  +  £A  (1  -  **)  +  e2m~2^W^  +  £2 


K^?)] 


+  Vi  [ee  (  U,p  +  +  er+m  f  ,PW,P  +  e2m  +  yy\  (7.23) 

,  -2*  ^  kp  +  £j  ,  (ge  ~  X7)2  4-ffe  +  F)2  ^  ]  j 

The  last  two  off-diagonal  constitutive  relations  for  the  out-of-plane  stress  components  yield  the  scaled  rela¬ 
tions: 

e^iSnzi  =  y  [e^1 1?*  +  e”W,p  +  er+’n~1f tPWx  +  e^W,^  +  e2^1  ( U,PUX  +  VpVl( )] ,  (7.24) 


Ve  +  u\  ,  _r+m  ft  ^  .  ,2m  / 1  f??2  .  ^,1 


+  er+mrpPTp  +  e2m  ijy2„  + 


epSiSezi  =  $  |e'_1V;c  +  £m  —  +  ^2m  1  *4“  £  f - - 


?  (7.25) 


respectively,  and  the  scaled  version  of  the  off-diagonal  in-plane  constitutive  relation  (5.6)  has  the  form 

=  Si  h  U  +  S*zl\  „  LMsl  .  EeE* 


m  Fp(t/,e-y)  +  yp(Fe  +  g) 
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The  final  two  constitutive  relations  (5.3)  and  (5.4)  take  the  following  forms: 
en^iSRRi  =  enXiS?m  +  Si  f  (1  -  Vi)  { elUtP  +  er+mT>pW,p  +  \  [e2mW2p  +  e21  (ft2  +  V2)] } 

+  em~xW^  +  e2m~2  ~W2(  +  e2t~2  ^  (ft2c  +  V; 


(fte-VQ2  +  (V(e  +  ft)2 
2  P2 


■  (7.27) 


and 


w2 

,  r2m  .e  ,  2 1 

2  p2  + 


snHiSeei  =  enXiS?m  +  £i  ^(1  -  Vi)  je* 

+  U,p  +  er+mT,pWp  +  \  [e2mW2  +  e21  (ft2  +  V2  )] 

I(^<)})- 


(ft,e-VQ2  +  (fr,e  +  ft)2 

2p2 


11 


+  em~vWtQ  +  e2m~2 -W2(  +  e2l~2 


(7.28) 


Next,  applying  the  scalings  obtained  in  (7.1)— (7.5),  (7.6)-(7.8),  and  (7.9)-(7.14)  to  the  right-hand  sides 
of  equations  (4.16)-(4.24),  we  obtain  the  stress  components  Tabi  in  terms  of  our  scaled  variables: 


Trri  —  en  EiSRRi  +  en+e  Ej 
Tr&{  =  en  SfSi?©i  +  en+*  Hi 

TRZi  =  epXi  SRZi  +  e^Si 


i  (u,pSRRi  +  Sfl©^  +  et+p~1  Si  U,cSRZi, 

!»  ^,P5fl©i  +  ^’e~  y  5eei^  +  el+p~l  S iUxSezi, 

i  (u,pSRZi  +  U'e  ~  V Sqz^J  +  e/+9_1  E iUjSzzt, 

Tern  =  enXiSRei  +  en+l^  (vJRRl  +  ^jEsRe^J  +et+p~1  E iVxSRZi, 

Teei  =  enXiSeei  +  £n+t^i  (vjnei  +  ^y^-See^J  +et+p~1  E iVt(Sezi, 

TeZi  =  epXiSeZi  +  e^Si  (v,pSRZi  +  - 5eZij  +e<+9"1  XiVi(SZZi, 

W,B  ; 


(7.29) 

(7.30) 

(7.31) 

(7.32) 

(7.33) 

(7.34) 


=  ep T.iSnzi  +  er+n  E;  fpSRRi  +  em+n  Si  |WlPSfl*i  +  +  em+p_1  S^S^i,  (7.35) 

Tz©i  =  ep  S iSezi  +  £r+n  S<  r,pSftei  +  em+n  Ei  fepSfl©i  +  ^SeeA  +  em+p~l  E iWxSeZi,  (7.36) 


T^zi  =  e«  S iSzzi  +  £r+p Si f ,pSRZi  +  em+p  Si  ( W<pSRZi  +  ^Sezi )  +  em+fl-1  SiW,c5^i.  (7.37) 

V  P 
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Substituting  these  scaled  expressions  in  the  equilibrium  equations  (4,12)-(4.14)  yields  equilibrium  equations 
in  terms  of  the  scaled  components  of  the  second  Piola-Kirchhoff  stress  tensor: 


en  SRRi,p  +  en+l  I  UJRRi  +  U'9  +  e 


1  (ut<sRZi) 


+ 


Ue-V: 


£n  SjiQi^Q  +  £n+ *  UtPSR@i  H  —  ~5©©i  +  £ 


fft+p-i 


+  ep  1  Snzi£  +  £P+l  1  UiPSRZi-\ — — - Sezi  +  £J 


}© 


U^e-V*  \  , 


(Ujezi)€ 

(Ujzzi) 


+  -  { en  SRm  +  en+i  ( uJRRi  +  ^’e  -Sfieil  +et+p~1  UxSRZi 


enSe&i  +  en+e  (vpSRet  +  ^e^5eeJ  +e<+p"1  VASezi 


=  0, 


(738) 


enSR@itP  +  en+i  VpSRRi  +  V'e  +  USR@i j  +£i+p  1  ^cSrzi) 


+ 


en  Seew  +  en+e  (vjaei  +  +  e/+p_1  (^<Sezi)  Q 


^ ^+?-2  (y^Szzi) 


+  e*"1  Sezi,(  +  +  r’*>J"SBziJ  +  et+9~ 2 

+  ^ | en  5sei  +  £n+*  (VpSRRt  +  -SRei)  +£t+p~1VxSRzi 


X 


+ 


sn  SRQi  +  £n+i  (uJRBi  +  -See*  J.  +  et+p  1  U^Sezi 


=  05 


+  em+n 


SRRi  +  ^SRei )  +  £m+p_1  (wxSRZi) 

pi  \  /  *P 

/  j  P 

WBl 


+ 


(739) 


epSRzi,P  +  er+n  (f  ,pSjiih) 

5 Sezi,s  +  £r+n  (f,PSRei)  e  +  em+n  {wJrsi  +  ^5©0ij  +  £m+p“1  (^ezi)  ( 

+  e9"1  Szzix  +  er+p~1  (f  JRZi)  c  +  e!n+p’1  (wJRZi  +  ^Sez*]  +  £m+9“2  (wjzzi)  ( 


epSRZi  +  £r+n  r,PsRRi  +  £m+n  (  W,pSRRi  +  ^SHei  ]  +  em+p~l  WASRZi 


+  =  0, 


where  the  gravitational  body  force  has  been  assumed,  as  stated  earlier,  to  scale  as  follows: 

apotg  =  2i£*5- 


(7.40) 

(7.41) 
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Finally,  assuming  that  the  pressures  p*  scale  like  the  stress  tensor  component  Szz,  i.e., 


- 


(7.42) 


our  boundary  conditions  (6.5)-(6.7)  take  the  following  scaled  forms,  after  eliminating  common  factors  of  E* 
on  both  sides  of  each  equation: 


(7.43) 


epS£z  +  e»l 


V>P 


Si z  +  s|2  +  v±sfz 

gjgg-jfefg; j-fo) 


=  e«+mp±— ^  +  e«+m+* 


p*. 


(7.44) 


£"  Sfz  +  £m+'  I  W'iS*7  +  )  +  £q+m_1  H'JSf 


e-9+2^ 


+  u* 


u* 


(7.45) 


The  method  of  asymptotic  expansions  proceeds  from  this  point  by  introducing  asymptotic  series  expan¬ 
sions  in  a  new  parameter  8  =  e**  (where  \i  is  yet  another  exponent)  for  each  of  the  scaled  stress  tensor 
components,  as  well  as  each  of  the  scaled  displacement  components  and  its  partial  derivatives.  Thus,  denot¬ 
ing  by  x  and  y  any  two  such  components  (or  partial  derivatives  of  displacement  components),  we  set 


x  =  ^2skx{k),  y  =  J2skV(k)>  (<5  =  £m). 

k= 0  k=0 


(7.46) 


Various  theories  of  two-dimensional-like  bodies  are  obtained  by  making  particular,  essentially  ad  hoc ,  choices 
of  the  scaling  exponents  in  equations  (7.2)~(7.5)  and  equation  (7.46).  For  example,  if  we  follow  Tarn  [12]  and 
set  r  =  m  =  1,  £  =  n  =  2,  t  =  3,  p  =  3,  q  =  4,  and  \i  -  2,  we  obtain  a  generalized  geometrically  nonlinear, 
laminate  shell  theory  (equivalent  to  von  K&rm&n  plate  theory  when  specialized  to  a  single  material  that  is 
initially  fiat,  so  that  hc  =  0  and  T(R)  =  0,  and  gravitational  body  forces  are  ignored).  If,  on  the  other  hand, 
we  follow  the  work  of  Erbay  [11],  and  set  r  =  m  —  1/2,  l  =  n  =  1,  t  =  3/2,  p  —  2,  q  =  5/2,  and  p  =  1/2, 
the  leading  order  equations  are  a  generalization  of  geometrically  nonlinear  membrane  theory  (equivalent  to 
Hencky-Campbell  membrane  theory  [13,  14]  when  specialized  to  a  single  material  that  is  initially  flat,  and 
gravitational  body  forces  are  ignored).  We  have  also  found  that  by  setting  r  =  1,  m  =  2,  £  =  n  =  3,  t  =  4, 
p  =  4,  q  =  5}  and  p  =  1,  the  leading  order  results  are  those  of  a  geometrically  linear,  laminate  shell  theory 
(equivalent  to  classical  laminate  theory  [15]  when  F(JF2)  =  0),  while  r  =  m  =  l,  £  =  n  =  3/2,  t  =  5/2,  p  =  3, 
q  =  7/2,  and  p  =  1/2,  yields  a  theory  of  a  geometrically  linear  membrane  laminate.  These  characterizations 
of  various  theories  by  particular  choices  of  the  scaling  exponents  are  tabulated  in  Table  1.  In  the  Sections 
that  follow,  we  develop  each  theory  based  on  the  choice  of  exponents  indicated  in  this  Table,  where  the  terms 
” shell”  and  ’’membrane”  distinguish  the  dependence  of  the  leading  order  in-plane  displacement  components 
on  the  axial  coordinate  Z :  for  a  membrane  they  are  independent  of  Z,  for  a  shell  they  axe  linear  in  Z. 
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Table  1;  Values  of  sealing  exponents  for  various  theories. 


8  Geometrically  Nonlinear  Shell  Laminate  Theory 


Following  the  work  of  Tarn  [12]  we  set 

r  =  m  =  1, 

£  ~  n  =  2, 

£  =  p  =  3,  and  g  =  4 

in  equations  (7.5 

>)-(7.5). 

Thus,  we  have 

r  =  eaT,  Uz 

=  saW , 

(7k  =  £2oE7,  Dre  =  £2aV, 

(8.1) 

Sum  = 

E  S%Qi  — 

Suet.  = 

=  £2Si5fl0i,  ST 

=  e2SjSfm, 

(8.2) 

Srz%  - 

s^iSuzu 

S&zi  = 

£3SiS0£j, 

(8.3) 

Szzi  = 

U  P±  = 

e4S±p±, 

opoi  S  =  Si  £3  5- 

(8.4) 

8.1  Leading  Order  Results  Obtained  by  Scaling  of  the  Constitutive  Relations 

Beginning  with  the  constitutive  relation  (7.23)  for  Szz,  we  obtain 


e^iSzzi  =  +  Ed  (1  -  Vi) 


w*  +  \k +e2  ( ^4-^ 


+  ^[e2  ft;,p  +  f,pW,p  +  iw2p  +  ^^  +  0 


+  0(e4). 


(8.5) 


Again  following  Tarn  [12],  we  set  p  =  2  in  equation  (7.46),  and  find  for  the  product  of  any  two  of  the 
asymptotic  expansions  (7.46),  to  tenth  order  in  e  (recalling  that  5  =  =  s2): 

*5  =  S  a(*i+fca)  *(*>£(*,), 

k\=  0  fc2=0 

=  *(0)5(0)  +  [*(o)5(i)  +  *(i)5(o)]  e2  +  [*(0)5(2)  +  *(i)5(i)  +  *(2)5(0)]  e4 

r^.^.  ^  ^  ^  ^  (8.6) 

+  l*(0)2/(3)  +  *(1)»(2)  +  *(2)P(1)  +  *(3)J/(0)J  £ 

+  [*(0)5(4)  +  *(1)5(3)  +  *(2)5(2)  +  *(3)5(1)  +  *(4)5(0)]  e% 

+  [*(0)5(5)  +  *(1)5(4)  +  *(2)5(3)  +  *(3)5(2)  +  *(4)5(1)  +  %)5(0)]  £10  +  0(£12). 
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Substituting  in  equation  (8.5)  the  asymptotic  series  for  each  of  the  variables,  we  obtain 

s'E.Smzz,  =  S!s,sr  +  £,j(l- n)  Wm,(  +  \w}ou  +  C  f  Hf(i).c  +  %),< 


t 


-)]+-Ke« 


(0  ),p  +  F,PW(0),p  +  2W(0),p  (8-7) 


V(o),e  +  U(Q)  % 


+ 


g5)] }  +  <V)- 


P  2  p* 

The  leading  order  term  of  this  relation  yields  our  first  important  result: 

tiW  +  (  s  \%u  (2  +  c)  =  0. 


(8.8) 


Thus,  W(0),c  must  satisfy  either  W(o),f  =  0,  or  W(o),c  =  -2.  In  order  to  eliminate  the  second  possibility,  we 
appeal  to  the  form  of  the  Jacobian  determinant  under  the  scalings  being  considered,  viz.,  from  (7.22): 


J  =  l  +  wA  +  e2 


-  Ve  +  U 
U,p  +  -----  +  W. 
P 


;,C  (u,p 


+ 


V}e±U 


v<  -  r,pi/,c 


+  £ 


u 


,p 


V#±u 

p 


B-V 


v^_±u_ 

p 


V\e-V 


)-'-P 

(v,pwx  -  wtPvx  -  f,p  vx)  +  (viPu,<  -  V,p% c) 


)(W-W?A-f,pf/,c)  (8-9) 


In  the  limit  e  ->  0,  the  condition  J  >  0  implies  that  we  must  have  W(0),£  >  -1,  which  precludes  the  second 
possible  solution  W(o),c  =  —2.  Thus,  we  must  have 


W(o),c  =  0  =*>  W(0)  =  w(p,Q), 


(8.10) 


where  w  is  an  arbitrary  function  of  p  and  0  only.  Under  these  conditions,  equation  (8.7)  reduces,  after 
dividing  through  by  e2,  to 


e2S«5(  0)zzi  =  XiS?m  +  £i 


(1  -  Vi)  W(1),c  + 


Vi  (u(0), 


‘) 


p  +  t,p™,p  +  2™.p  + 


1-2  ^w.e  +  U(0) 

P  2  p2 


+  0(e 2) 


(8.11) 


The  leading  order  term  on  the  right-hand  side  of  this  equation  yields  another  relation  that  will  be  needed 
later: 


W(1),C  + 


£(0 U±goU  _  XjST  « 


1-2  ,  1^(0)  ,e  +  ^(0) 


tu2e\ 

#J  •  <8-12> 


Next,  consider  the  two  off-diagonal  constitutive  relations  (7.24)  and  (7.25)  for  the  out-of-plane  stress 
components.  With  the  present  scalings,  we  obtain 

e^iSnzi  =  f  [(£,<  +  WtP  +  T<PWX  +  WjW<  )  +  £2  (  U,PUX  +  V,PVA  )] ,  (8.13) 
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(8.14) 


e2T,iSezi 


Gi 

2 


p  t  W,e  |  W^W^e 

,C  P  P 


+  e 


Uc(U,6-V)  +  VdVe  +  U) 


where  a  common  factor  of_e  was  cancelled  in  both  expressions.  From  the  leading  order  terms  of  these 
expressions,  recalling  that  W(o),c  =  0,  we  obtain  the  following  two  equations: 

•"v  —  U)  A 

U(o)x  +  ™,p  =  0,  V(o)X  +  =  0, 

which  can  be  integrated  to  obtain 

U(o)  =  2(p,0)  -  Cw,P,  and  F(0)  =  v(p,Q)  -  C^y ,  (8.15) 

where  u  and  v  are  arbitrary  functions  of  p  and  0  only.  Thus,  the  leading  order  results  have  provided 
non-dimensional  forms  of  the  well  known  Kirchhoff-Love  expressions  for  the  displacement  components. 

The  scaled  version  of  the  off-diagonal  in-plane  constitutive  relation  (7.26)  reduces  under  the  present 
scalings  to 


Gi 


V  f  k  LeE*  t  e2  up(Ue-y)  +  vAye  +  u) 

,p  p  p  p  *  * 


where  a  common  factor  of  e2  was  cancelled.  Thus,  to  leading  order  we  obtain 

^  c  _  Gi  (tr  ,  ^(o),e  “  V(o)  ,  f,pwe  w,pw,e ' 

ZiS(o)Rei  =  y  1 vm>p  + - - - +  — —  +  -syL- 

Introducing  an  in-plane  strain  component  crq  defined  by 

-  -  1  (y  ,  G(0),e  -  V{0)  rp«5,e  $Pw,e\ 

eRe ~ 2 [VW'<>+  p —  ~r~)' 


„  |  g,e-g  [  r^w.e  t  w,pw,Q  ( W,pe  w,e\ 

P  P  P  \  P  P2  ) 


=  £j?e  -  {hue, 

where  we  have  introduced  (-independent,  terms  in  the  last  line  defined  by 


^0  _  1  ,  ue-v  r,pw,0  ,  w,pw,e\ 

£fle  =  2  (  +  ~T~  ~P~  ~P~ ) 

we  can  write  this  constitutive  relation  as  simply 

2 iS(G)RBi  =  GitRQ. 


T  —  w,pB  w  © 

™R&  -  o  ? 

p  p 2 


(8.16) 


(8.17) 


(8.18) 


The  final  two  constitutive  relations  (7.27)  and  (7.28)  take  the  following  forms  under  the  present  scalings: 


e2XiSRRi  =  e2XiS?m  +  £  j  (1  -  Vi)  e2  (  UiP  +  f ,PW,P  +  \w2p  ) 


+  Vi 


^  w&\  ^  ^  ,  ifj>2 ,  j  r 


+  I  +  Wx  +  2W *  +  £' 


2p 


(8.19) 


+  0(eA), 
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and 


e2S iSeei  =  e2Si5fm  +  SU  (1  -  Vi)e 


„2  / ge±_g  ^  W?e 


+  Je2  I  U,p  +  f,pW,p  +  \wl )+Wx  +  ^  +  e2 


+  0(e4). 


The  leading  order  terms  involving  W ^  in  these  expressions  vanish  according  to  (8.8),  and  (8.12)  can  be  used 
to  replace  the  second-order  terms  involving  W^(i),{,  yielding 

ys.  ^  f  ^  ,  Q71T71  jyr.  I  I 

ZiS(o)RRi  =  ViS?™  +  £il  (1  -  i 'i)eRR  +  Vi  ?ee  -  ^  ^rr  +  ?e0)  r>  (8-21) 


(l -Vi) tee  +  Vi  eRR  -  ^  z ~t i %  ~  T^~Vi  ^ RR  + 
where  we  have  introduced  scaled  in-plane  strain  components  ?rr  and  e©©  defined  by 

eRR  =  U(o),P  +  f  ,pw,p  +  \w2p  =  uiP  +  f iPw,p  +  \w2„  -  (w,pp  =  e°RR-  (kRR, 


__  V(o),e  +  U(o)  w%  __  vye  +  ^  ,  w>e  r  ( ^ 


Cee  -  n  +  9n2  n  +  9^2  ^  V  n  + 


^ee\  _ 

P2 )  - 


e@©  —  (fcee- 


The  (-independent  terms  appearing  in  the  last  two  equations  are  defined  by 

?RR  =  «,p  +  f ,pWp  +  \<V2p  ,  kRR  =  W,pp, 


<o  _  v,e+u 
eee  =  : 


_  w,p  ,  u>,ee 


+  7TT  >  ^ee  =  —7-  + 


Replacing  £<  by  its  definition  (5.7)  in  (8.21)  and  (8.22),  and  simplifying  the  results,  yields 

^  ^  E' 

E iS(0)RRi  ~  d-  ^  _  ^2  (  ^RR  ^te0©  )  > 


^i^O)©©*  =  ^>iSi  +  _*  2  (?©©  +  VitRR)  , 


where 


a  =  i^§r- 

1  -  Vi 


We  conclude  this  Section  by  rewriting  the  important  results  in  terms  of  leading  order  variables,  that  are 
functions  of  the  physical  coordinates  R ,  0,  and  Z,  viz., 


u  =  e2au ,  v  =  e2av,  w  =  saw 
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(8.31) 


U(o)Z  =  eaW{ 0)  =  w,  U{0jR  =  e2aU{0)  =  u  -  Zw,R,  U(0)e  =  e2aV{0)  =  v  -  Z^§-, 


R 


S(o)nei  =  £2EiS(o)nei  =  Giene, 


(8.32) 


S(0)RRi  =  £  ZiS(0)RRi  =  Si  +  -j~~2  (e-R-R  +  We©  ) 


(8.33) 


Ei 


%>)e©i  =  t  S,S(0)©©i  =  Sj+  -  — -  ?  (fee  +  v^rr )  ■ 

■*■  v% 


(8.34) 


where 


1 

eRe  =  e  Cfi©  =  - 


v'R+~ir+^r  + 


m m 

w,®\ 

R 

\  R 

R?  ) 

4e  -  (8.35) 


=  £  crr  =  uR  +  TiRwtR  +  -vfR  -  Z  wtRR  =  e°RR  -  Z  kRR, 


(8.36) 


_  2-  V,©  +  U 

eee  =  eeee  =  ^  +  ~ 


,  w,ee\ 


+  “grj  =  4©  -  2*ee, 


(837) 


the  ^-independent  terms  of  the  last  three  equations  are  given  by 
1 


w,Rwt&\  _  w,Re  w.s 

+  R  )’  "A©  =  l~  ~ft2> 


(8.38) 


4a  =  u,r  +  r  ,rW'R  +  -w2r  ,  kRR  —  wtRR , 


(8.39) 


,o  _  w,©  + u  , 

00  “  R  +  2R2 


kee 


wji  ,  w,ee 
R  R2  ’ 


(8.40) 


and  in  (8.33)  and  (8.34)  we  have  introduced  the  in-plane  residual  stresses  defined  by 


Si  =  e2  S,-  Si 


1  —  2V{  ©  1  —  2z4- 

_ i^V.cnm  _  *  nnm  «.m 

i-*  A  -T3^si  - 


(8.41) 


where  the  last  equality  of  (8.41)  follows  from  (5.2),  and  refers  to  the  particular  case  where  the  residual 
stresses  are  purely  thermal  in  origin. 

Note  that  only  three  in-plane  constitutive  equations  (8.32),  (8.33),  and  (8.34)  are  obtained  as  leading 
order  results,  from  which  the  associated  in-plane  stress  components  can  be  determined  in  terms  of  the 
displacement  components.  The  other  constitutive  relations  provide  the  leading  order  (Kirchhoff-Love)  forms 
(8.31)  of  the  displacement  components,  but  the  associated  out-of-plane  stress  components  will  be  shown 
later  to  be  determined  in  terms  of  the  in-plane  stresses  via  first  integrals  of  the  equilibrium  equations. 
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8.2  Equilibrium  Equations  to  Leading  Order 

Under  the  scalings  of  this  Section,  the  equilibrium  equations  (7.38)-(7.40)  can  be  written  as 
£2  Snm,p  4-  e4  (U'p§RRi  +  — : +  e4  (0^SRzi'j 


) 


,p 


1 

+  — 
9 


+  £4  f  U,pSRei  4 — ^ — S'©©!]  +  £4 


(u,p 


,© 


,© 


U,e-V~ 


4  e2  SRZix  +  e4  1  +  —^—Sezi  J  +  £4  (£,cSzzi)  ( 


.C 


P  l  V  P 


+  *(e2&uii  4-  e4  (uJRRi  +  U£-V-SRe?J  +e4UxSRZi 


£2  Seei  +  £4  ^pSflei  4 — — S©©j^  4-  £4  V^Sezi  ^  =  0, 


(8.42) 


£2  Sjieitp  4*  £4  (  + 


U,e  + 


-SRei)  +£4  (vxSRzi)f 

'  >P  i 


1 

4-  - 
9 


£2  S&&i, e  4- £4  ^UiPSfi©j  4-  ,e^ — 5©©i^  4-e4  (v^Sezi^J 


,© 


.© 


+  e2  Sezi£  +  e4  (^PS*zi  +  ^±-^Sezi  J  +  e4  (v>cSZZi)  ^ 

i|e25flei  +  £4  (v,pSRRi  +  +  £4  U<SflZi 

(u,pSRei  +  ^e~^,See,j  +  e4  U,cS©zi  j  =  0, 


£  SRGi  +  £  (  U,pSRei  + 


(8.43) 


£3  SRzi,p  +  £3  piSflfli  j  ^  +  £3  H ^-SRQiJ  +£3  ^W^S/jz;  j 


.p 


1 

+  - 
p 


+  e1 


1 

+  - 
P 


£3  Sezi,e  +  £3  ^r,p5fl©i j  ^  4-  e3  ^VU)P5fl©j  4 — “Seei  ]  +  £ 


*) 


1  (w,<Sez<) 


,© 


.© 


Szzi.c  +  £3  (f  ,pSRZi)  c  +  £3  f  WJRZi  +  -  fSezi  J  +  £3  (w^Szzi) 

'  '  iC 

1  SRZi  +  £3  ftPSRRi  +  e3  +  ^S*©ij  4-  £3  WxSRZi 


< 

+  £3  9  =  0. 


To  leading  order,  using  the  earlier  result  W(o),(  =  =  0  in  the  third  equation,  we  thus  obtain 

S(.0)RRi,p  H - S(o)fl©i,e  4-  S(o)flZ«,<  H - (S(o)flflt  —  S(0)©©j j  =  0, 


(8.44) 


(8.45) 
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(8.46) 


S(o  )RQi,p  +  ^5(o)ee»,e  +  Scoezi.c  +  -S(o)nm  =  0, 
u,pS(o)Rm  +  -f-S(0)Rei  4-  S(o)fi2i]  p  +  ~n  wiPS( Q)Rei  +  —5(0)00,.  +  5(o) ezi 

H  p  p  J,e 

+  U!pS(o)RZi  +  —  S(0)@Zi  +  S(Q)ZZI  +-  ^,pS(0)RRi  +  +  S^0)RZi  +5  =  0, 

where  we  have  introduced  in  the  last  equation  a  new  p  and  ©-dependent  function  S3  defined  by 

w(p,0)  =  w(p,  @)  +  f(p). 


8.3  Scaled  Boundary  Conditions  of  Pressure 

The  boundary  conditions  of  pressure,  equations  (7.43)— (7.45),  are  given  under  the  present  scalings  by 


e3S±z+s5  USSiz+f-*  ^ 


Siz  +e5U±S % 


=  e5p±W f 


7  ty*  (fg + e*y 


e35±  +e5  F±5±  + 


vi  +  u± 


Siz  +  c5t>fS|z 


=  ,»s±Et  ^  ,7  [  ^-^(g.e-^) 


£45|z  +  e4  ( l+±5±  +  -^5±  ]  +  ^mst 


-  -zz  -r  ^  ^",p+rz  T  —r^ezj  f  e  ^,c^zz 

=  -'*?*-  e“  (vi  +  p*  -  <■  ej 

To  leading  order,  they  reduce  to 


v±  +  u± 


.9*  (GAzll 
l  p± 


(8.50) 


c± 

°(0  )RZ 


c± 

°(0)©Z 


ci  _  £>± 

°{0)^Z  “  ~P  * 


where  the  derivation  of  the  last  result  for  S^Zz  squired  the  previous  two  results  for  S^RZ  and  S^e^, 

and  the  fact  that  ^q)»c  =  Stated  in  terms  of  physical  functions  of  the  physical  coordinates,  these  take 
the  forms 


5(0)JJZ  -  °>  5(o)ez  -  °>  5(0)zz  =  e4s±%)zz  =  =  -p±.  (8.53) 

8.4  Leading  Order  Equilibrium  Solutions  for  the  Out-of-Plane  Stress  Components 

In  terms  of  leading  order  physical  functions  of  the  physical  coordinates  R,  0,  and  Z,  the  leading  order 
equilibrium  equations  (8.45)-(8.47)  are  given  by: 

S{0)RZi,Z  +  ‘S'(0)Rfli,«  +  -jjj  {S(o)RQi,e  +  S(0)RRi  ~  <S(0)©©,)  =  0,  (8.54) 
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5(o)ezi,z  +  S(o)Rei,R  +  ^  (5(o)©ei,e  +  2  5(0)r©i)  —  0, 

+  -j£-S(0)ezi  +  -5(o)zz<]  z  +  \  [w,fl5(o)iie<  +  -^^(o)©©;  +  5(o)ez»]  € 

1  *  -i 


+ 

where 


XI  J  at  L  At,  - 

[w, RS(0)RRi  +  -j^-S(0)nei  +  5(0)flZi]ifl  +  ^  [<*>, RS(0)RRi  +  -^5(0)r© i  +  £(0)/^]  +  poi  5 


w(R,0)  =eau{p,Q)  =  w{R,&)  4-  T(i?). 
Equations  (8.54)-(8.56)  can  be  rewritten  as 

S(o)RZi,z  +  ^  ^(.RS(o)flRt)  R  —  5(o) ©©i  +  5(o)R©i,©j  =  0, 


1  /^2  -  '  * 


5(o)©z«,z  +  -^2  (-R2  5(0)R©t),fl  +  —S( o)©©<,©  -  0. 


(8.55) 

(8.56) 

(8.57) 

(8.58) 

(8.59) 


and 


[utRS(Q)Rzi+-j£-S(.o)ezi  +  5(o)zz<]  z  +  R  [w.-R^(o)«©«  +  ~^5(0)©©i  +  5(0)©z>]  Q 

+  ^  (H5(o)rz<)iR  +  ^  [-R  (w,R5(0)RRi  +  -jpS(0)fi©i)]  R  +  Poi9  =  0, 


(8.60) 


respectively.  Now,  according  to  (8.32)-(8.40),  S(o)j?©i?  £(o)itRi>  and  S(o)©©i  are  linear  functions  of  Z,  with 
coefficients  depending  only  on  R  and  0,  hence  equations  (8.58)  and  (8.59)  can  be  written  as 


S(o)RZi,z  +  oo*(-R>  ©)  +  Z  au(R,  0)  —  0,  S(o)&Ziyz  +  boi(R,  0)  +  Z  0)  —  0,  (8.61) 

respectively,  where  aoi,  an ,  boi,  and  bit  are  rather  complicated  Z-independent  functions.  These  two  equations 
can  be  easily  solved  to  obtain 

Z 2 

5(o  )RZi  +  Zaoi  +  —au  =  Frz«,  5(0)©zi  +  Zb0i  +  —  bu  =  Fezu  (8.62) 


where  Frz%  and  Fezi  are  arbitrary  functions  of  R  and  0  only.  Applying  the  first  two  boundary  conditions 
of  (8.53),  we  obtain  from  the  first  and  second  equations  of  (8.62): 

Frzs  =  —  ao«  +  "g*  o.is,  Frzc  =  —  g  °0c  +  "g"  °ic,  (8.63) 

^  „  h.  h\  _  h,  h2.  ,*&a\ 

F©z«  =  2b°s  +  ~8  FeZc  ~  ~  2  b°c  +  ~8  blc’  8-64 

respectively.  For  a  two-layer  laminate,  the  solutions  for  the  stress  components  in  the  two  materials  can  thus 
be  written  as 

S(0)«r,  =  (f  -  z)  00.  +  l  (*!  -  4ZJ)  «...  Sm„z,  =  -  (j  +  z)  «o,  +  j  (fc!  -  4Z!)  alc,  (8.65) 

=  (|  -  z)  k,  +  l  (d2  -  4Z!)  h„  Smbz,  =  -  (5  +  z)  Soc  4  i  {kJ  -  4Za)  bic,  (8.66) 
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and 

5(o)©z« 


respectively,  where  i  =  1  =  c  denotes  the  coating  and  i  =  2  =  s  denotes  the  membrane  substrate. 

Next,  we  note  that  since  S(o)HZi  and  5(o)ezi  are  quadratic  in  Z,  the  axial  equilibrium  equation  (8.60) 
can  be  written  as 

[w.RS(o)jizj  +  z  +  coi(R, 0)  +  Z cu{R.  0)  +  Z2c2i(R,Q)  =  0,  (8.67) 

where  coi,  ch,  and  c2;  are  complicated  functions  of  R  and  0  only  (note  that  Coi  includes  the  term  poi  g). 
The  general  solution  of  (8.67)  is 

w  ©  Z2  Z3 

u,nS(o)RZi  +  ~^-S(o)ezi  +  S(o)zzi  +  Zcoi(R,  0)  +  —  cu(R,  ©)  +  —  c2i{R,  ©)  =  Fzzi(R,Q),  (8.68) 

where  Fzzi  is  another  arbitrary  function  of  R  and  ©  only.  Applying  the  three  boundary  conditions  of  (8.53), 
we  obtain 

„  +  h  h 2  h 3  *  „  ■  -  h  h2  h3 

tzzs  =  +  -cqs  +  —  cls  +  —  c2s,  and  FZZc  =  ~P  -  ^  °oc  +  -g-elc  -  —  c2c.  (8.69) 

The  solutions  for  this  stress  component  are  thus 


s(o )zzs  =  ~P+  +  ~  Z'j  cos  +  g  {h2  -  4  Z2)  cls  +  ^  (h3  -  8  Z3)  c2s 


~  U,rS(q)rZs  ~ 


R  ‘ 


(8.70) 


5(o  )ezs 


in  the  membrane  substrate,  and 


5(0  )ZZc 


~  (|  +  z)  coc  +  |  ( h 2  -  4 Z2)  cle  -  ±  (h3  +  8 Z3)  c2c 


u,rS(o)rzc  -  ~^-S(0)ezc 


(8.71) 


in  the  coating. 


8.5  Continuity  Conditions  on  the  Out-of-Plane  Stress  Components 

Requiring  that  the  stress  components  S(0)RZ  and  S(0)ez6  be  continuous  at  the  interface  Zj  =  (hc  -  hs)/2, 
we  find  from  (8.65)  and  (8.66)  that  we  must  have 


and 


S(o)rzs(Zi)  -  S(0)RZc(Zi)  =  hsa0s  +  hca0c  +  ^hche  (als  -  oic)  =  0, 


(8.72) 


S(o)ezs(Zi )  ~  S(0)&Zc{Zi)  -  hsb0s  +  hcb0c  +  i hchs  (6is  -  bic )  =  0,  (8.73) 

respectively.  Assuming  these  continuity  conditions  to  hold,  the  continuity  condition  for  S(0)zz  at  the  interface 
follows  from  equations  (8.70)  and  (8.71),  viz., 


S(o)zzs(Zi)  ~  S(o)zzc(Zi)  =  p  +  hscos  +  hccoc  +  (ci*  -  Cic) 

l  1  (8.74) 

+  12  hi)  c2«  +  he  (3 h2s  +  hi)  c2c  =  0, 

where 

p  =  p~  -  p+  (8.75) 
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is  the  pressure  difference  between  the  lower  and  upper  faces  of  the  deformed  configuration.  Equations  (8.72)- 
(8.74)  represent  conditions  that  must  be  satisfied  by  the  fourteen  Z-independent  functions  ooi,  an,  &o»,  &i*> 
coi,  cu,  and  cu  (i  =  c  or  s),  in  order  for  the  out-of-plane  stress  components  to  be  continuous  across  the 
coating/membrane  interface. 

To  facilitate  the  development  of  the  continuity  conditions,  we  first  separate  each  of  the  in-plane  consti¬ 
tutive  relations  (8.32)-(8.34)  into  its  ^-dependent  and  Z-independent  parts: 

S{o)Rei  =  crRei  ~  ZriRQi ,  5(o)fijii  =  <?Ri  -  ZrjRi,  5(0)601  =  a&i  -  Zrfei,  (8.76) 

where 

0Rei  =  Gi  ,  t]RQi  =  Gjfcfl©,  (8.77) 

o'flt  =  St  +  Qi  (^rr  +  Vi c©e)  t  WRi  —  Qi  (&. RR  +  VikQ&) ,  (8.78) 


<TQi  =  Si  +  Qi  (e©e  +  v^rr)  ,  =  Qi  (k&&  +  VikRR) .  (8.79) 

In  (8.78)  and  (8.79)  we  have  introduced  material  parameters  Qi ,  defined  in  terms  of  the  moduli  and  Poisson’s 
ratios  by 


Qs  = 


i-vV 


Qc  = 


Ec 


I-vV 


(8.80) 


and  from  (8.38)-(8.40),  the  Z-independent  “strains”  and  “curvatures”  are  given  in  terms  of  the  leading  order 
displacement  components  by 


4e  =  g  yv’R  + 


R 


■  + 


R 


+ 


R  )' 


R 


(8.81) 


eRR  =  U,R  +  r, flW.fi  +  ■z’UTr  ,  kRR  =  W,RR, 


(8.82) 


o  _  v,e  + u  ,  w,©  _  wtR  w,ee 

Cee  “  R  +  2R2  ’  ee  "  R  + 


R? 


(8.83) 


The  development  of  the  continuity  equations  first  requires  identification  of  the  fourteen  Z-independent 
functions  aoi,au,boi,bu,coi,cu,  and  cu  (i  =  c  or  s).  Comparison  of  (8.62)  with  equations  (8.58)  and  (8.59) 
yields 

(8.84) 

(8.85) 


oo,  +  Zaii  =  —  (RS(o)RRi)  R  ~  5(0)eei  +  5(o)fi©t,ej , 


boi  +  Zbu  ~  (R?  S(o)Rei),R  +  —  5(0)©ei,©. 


R? 


R 


Similarly,  comparison  of  (8.67)  with  (8.60)  yields 

cot  -I-  Zcij  +  Z2c 2i  =  [w,fl5(0)B©j  +  ^^5(o)©©»  +  5(0)©zij  q 

+  (RS(0)RZi)  R  +  ^  [-R  (u,RS(o)RRi  +  -^•5(0)fi©i)j  r  +  poi9.  (8.86) 
Substituting  from  equations  (8.76)  in  the  right-hand  sides  of  (8.84)  and  (8.85),  we  easily  find 


aoi 


=  ^  [(■Rcrfl»),/i  _  °ei  +  dfl©i,©j  ,  an  =  [(7Z??fl»),fl  _  Wei  +  J?fl©i,©]  , 


(8.87) 
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bo»  -  jp  {R2  VRei),R  +  -^crei,©!  bit  =  —jp  (R2  r}Rei),R  - 


(8.88) 


for  either  i  =  c  or  s.  Substituting  from  equations  (8.76),  (8.65),  and  (8.66)  in  the  right-hand  side  of  (8.86), 
we  eventually  obtain  for  i  =  s: 


1  r  (  W'Q  h  h?  \ 

cos  =  Pos  g  +  [  I  w, RVRQs  +  +  2  hs  +  -g-Oi s  I 

/fin  h2  ta©  \  i 

+  ^2  RaOs  +  -gRals  +  RU'RCTRs  + -Jj-<TRQs  j  j, 

Cls  =  ~  +  “^^©s  +  bosj  e  +  (Eaos  +  Rj^,RT}Rs  +  -j^-TlRQs^j 


and  for  i  =  c: 


(8.89) 

(8.90) 

(8.91) 

(8.92) 

(8.93) 

(8.94) 

(8.95) 

(8.96) 

Beginning  with  the  continuity  condition  for  S(0)rz  given  in  (8.72),  we  substitute  from  (8.87)  into  that 
equation,  then  expand  the  derivatives  and  multiply  the  result  through  by  R,  to  obtain 


°2s  —  ~2jj  [bis,©  +  (-RaisXjjj  , 

-  If  f  _w  e  h ,  fi2 ,  \ 

Coc  —  POc  9  +  -jj  [  (  —  2  he  +  -g-bi e  1 


(  h7 

4 

[-21 

1  ‘ 

Clc  = 

R  . 

1 

C2c  — 

~2R 

8  J 


wte 
R  1 


hs  (R&RstR  4  GRs  —  &Bs  4  &RBste)  4  hc  {R&Rc,R  +  &Rc  ~  &Bc  4  <7R©e,©) 

—  2  hchs  [{RriRs^R  +  7]rs  —  7]Ss  4  7?H©s}©)  —  (R^Rc.R  4  *0Rc  ~  Vec  4  TfRBc.e)]  =  0. 

Substituting  from  equations  (8.77)-(8.79)  into  this  equation  yields,  assuming  the  residual  stresses  Sc  and  Ss 
to  be  constants, 


B(Ae°RRR  +  ^€©©,i?)  +  (A  —  Av)  ( e°RR  -  e©©)  +  Ae€%e  e 

+  R(BkRRiR  +  B^kee.R)  +  (B  —  Bu){kRR  —  fc©@)  4  -Be&i?©,©  =  0,  (8.97) 

where  we  have  introduced  the  following  constants  linear  and  quadratic,  respectively,  in  the  thicknesses  hc 
and  hs : 


A  —  hs  Qs  -j-  hcQcy  Av  —  ~f"  hgQci/£)  4©  —  hs  G s  *4“  h^Ga  (8.98) 

B  =  “  h€hs  {Qc  “  Qs )  ?  Bv  —  2  hchs  {Qc^c  ~  Qs^s)  Be  =  — hch8  (Gc  —  Gs).  (8.99) 

Similarly,  substituting  from  (8.88)  into  the  continuity  condition  (8.73)  for  £(0)©z,  we  obtain  after  some 
algebra  the  following  form  for  the  second  continuity  condition: 

Ae  (R2€RB)  r  4  i?(Ae©e  4  AueRR)  e  4  Be  ( R2kRe)R  4  R(Bkee  4  BvkRR)  B  =  0,  (8.100) 
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where  the  constants  are  defined  in  (8.98)  and  (8.99),  and  we  again  assumed  the  residual  stresses  to  be 
constants. 

The  final  continuity  condition  (8.74)  involves  some  rather  tedious  algebra,  which  we  omit  here.  However, 
in  order  to  write  the  result  in  terms  of  constants  previously  defined  in  (8.98)  and  (8.99),  we  remark  that  we 
make  use  of  the  following  identity: 

h\  Xos  ~  h2cxoc  =  (h*  —  hc)  (h*Xo«  +  hcxo c)  —  hch„  (xoc  —  Xo«)  >  (8.101) 

where  the  Z-independent  function  x  may  be  either  a  or  b.  Using  this  identity,  the  final  continuity  condition 
can  be  brought  to  the  form 

(he  _  M  [R  (-A  t°RR,R  +  ^eee,fl)  +  iA  ~  Av)  —  eee)  +  -Ae  £«©,©] 

—  [R  (B  e^R  R  +  J9„e©©  rt)  +  {B  —  Bv)  (crr  —  c©©)  +  Be  ei?e,e] 

-  [R  ( DkRR,R  +  DvkeetR)  +  (D  -  Dv)  ( kRR  -  fc©e)  +  DekRe,e] 

+  R  (Ar  +  A  crr  +  Av  e©©  +  B  kRR  +  Bv  fcee)  +  ^  (A©  eR ©  +  Bq  kRe)  ]  } 

+  -^2  ( h ,  -  hc)  ^4©  (R2eRe)R  +  ii(Ae©e  ©  +  A„e5efli©)j  -  [-B©  (R2e°Re)R 

+  R  (-Be©©,©  +  Bv  eRR  B)  j 

-  [l?©  (R2kRe)  R  +  R(Dkee,e  +  D„kRRie)^  +  #2  (-4©  e$i©  +  BekRe ) 

+  RW'G  (A/-  +  Ae©©  +  Av^rr  +  B  kee  +  BvkRR )  |  ^  +  R{p  +  7ofl)  =  0,  (8.102) 

where 

7o  =  hspos  +  hcpoc,  (8.103) 

is  the  areal  density  of  the  coated  membrane  (mass  per  unit  area,  Kg/m?,  of  the  circular  disk  perpendicular 
to  the  axis),  and  we  have  introduced  three  new  constants,  each  cubic  in  the  thickness  hs,  defined  by 

D  =  5|[(i  +  Zn)Qs  +  (3  +  n)H2Qc),  (8-104) 

XZ 

P„  =  £|[(1  +  3  %)Q.v.  +  (3  +  H)H2QCV c],  (8.105) 

XZ 

Be  =  J|[(l  +  3K)G*  +  (3  +  H)H2GC],  (8.106) 

as  well  as  a  constant  M  having  units  of  N/m: 

A r=hsSs  +  hcSc.  (8.107) 

In  the  definitions  (8.104)-(8.106)  we  have  eliminated  hc,  introducing  instead  the  thickness  ratio  H  defined 

by 

71  =  ^.  (8-108) 

h8 

The  terms  in  (8.102)  involving  the  constants  A,  A„,  and  Ae  can  be  eliminated  in  favor  of  terms  involving 
B,  Bv,  and  B@  and  three  new  constants  D,  P„,  and  De,  by  using  the  first  two  continuity  conditions  (8.97) 
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and  (8,100).  With  these  manipulations  the  third  continuity  equation  simplifies  a  bit  to 

{  -  [B(B^rr,r  +  +  (B  -  B„)  ( e°RR  -  e©e)  +  Bee?*©©] 

-  [R  (DkiiR'R  +  D„kee,R)  +  (B  —  Dv)(kRR  -  A;©©)  +  Befc#©,©] 

+  R  J  &,R  (A  4-  AcJjjj  +  Av  €@q  +  B  k-RR  +  Bv  A©©)  +  (A©  e°Rq  +  Bq  fcj*©)  j  J- 

+  ^2  {  ~  {^e  (R2£°Re)  'R  +  R  (-Be©©  ©  +  Bv€RRq) 

-  |b©  (R2kRe)  R  +  R(Dkee,e  +  Dv  fcnH,©)j  +  R2  utR  (A©  eRe  +  B©  kR&) 

+  Rwte  (AT  +  Ae©e  +  Ave°RR  +  BA:©©  +  BvkRR.)  j  +  R(p  +  'yod)  =  0,  (8.109) 

J  ,0 

where  D,  and  27©  are  defined  by 

D  =  ^{(1  +  M2)Qs  +  (3  +  H2)nQc],  (8.110) 

Du  =  ^[(i  +  3H2)Qsvs  +  (3  +  H2)HQcvc],  (8.111) 

Be  =  §[(1  +  3n2)Ga  +  (3  +  U2)%GC).  (8.112) 

We  conclude  this  Section  by  summarizing  the  continuity  conditions  in  terms  of  the  strains  and  curvatures 
(and  their  partial  derivatives): 

B  (A  £rr,r  +  Av  £©©>fi)  +  (A  -  A„)  (e°RR  -  e©e)  +  A©  eRe  © 

+  R{BkRRiR  +  Bykee^)  +  (B  -  Bv){kRR  -  fc©©)  +  SqA;#©,©  =  0,  (8.113) 

A©  (R2e° ie)R  +  R  (Ae©0  +  A„  cRR)e  +  B©  (R2kRe)  R  +  R(Bkee  +  BvkRR)  e  =  0,  (8.114) 

{  -  [B  (B  eRRR  +  Bve%e  R)  +  (B  -  B„)  (e°RR  -  e©0)  +  B©e^e  e] 

—  [R  (D  kRR<R  4-  Dv  &©©,i i)  +  ( D  —  Dv)  (kRR  —  fc©©)  4-  B©  fcj?©,©] 

+  B  [  (W,R  +  r,ii)  (Af  +  A  eRR  4-  Av  4e  +  B  kRR  4-  Bv  fc©©)  4 — (A©  e°R&  +  B©  fc/j©)  j  J 

+  ^2  {  _  [Be  (R2£Re)tR  +  R  (Be©e,e  +  BueRRe)j 
-  B©  ( R2kRe )  R  +  R  ( Bfc©©,©  +  B„  fciiR,©)J  +  B2  (tn,^  +  r^)  (A©  4©  +  B©  fcj?©) 

+  Uffi,e  {Af  +  Ae%e  +  Ave°RR  +  Bfc©@  +  BvkRR)  }  +  R(p  4- 7off)  =  0,  (8.115) 

where  the  definitions  of  the  constants  are  repeated  here,  for  convenience: 

M  =  h3Ss  4  hcSc ,  (8.116) 

-4  =  hsQs  4-  hcQ€l  Av  —  hsQsvs  4  hcQcuCj  A©  =  hsGs  4  hcGCj  (8.117) 

=  2^c^s  c  —  Qs)j  —  —heh8  (Qevc  —  Qsvs),  5©  =  —  (Ge  —  G*) ,  (8.118) 
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(8.119) 


0  =  M[(i  +  3  H2)Qa  +  (3  +  n^HQc], 

Dv  =  ^f[(l  +  m2)Q,vs  +  (3  +  n2)HQci'c],  <8-120) 

De  =  5|[(1  +  3?f2)G»  +  (3  +  H2)^Gc].  (8-121) 

Note  that  in  (8.115),  we  have  replaced  two  occurrences  of  the  partial  derivative  w,r  =  w.fl  +  F,r.  For 
comparisons  with  Wittrick’s  work  [16],  we  remark  that  our  w  corresponds  to  the  function  he  denotes  by  w, 
our  T  corresponds  to  the  function  he  denotes  by  w0,  and  our  w  corresponds  to  the  function  he  denotes i  by  w  . 
The  constants  given  above  differ  from  those  in  equations  (8)  of  [16]  due  to  a  difference  in  the  choice  of  origin 
of  coordinates  (our  origin  is  at  the  center  of  the  middle  plane  of  the  reference  placement,  while  Wittrick  s 

corresponds  to  the  center  of  the  interface  plane  between  coating  and  substrate  of  the  reference  placement). 

The  constants  Aq,  Be,  and  De  do  not  appear  in  Wittrick’s  paper,  since  he  treats  only  the  arisymmetnc 

problem. 


8.6  Formulation  of  Equilibrium  Equations  in  Terms  of  Stress  Resultants  and  Stress  Couples 

It  is  more  common  to  find  the  equilibrium  equations  for  a  shell  or  plate  presented  in  terms  of  stress  resultants 
and  couples,  and  in  many  respects  such  a  formulation  is  simpler  than  (though  equivalent  to)  the  one  discussed 
in  the  previous  two  Subsections.  However,  it  was  felt  that  the  approach  via  continuity  conditions  was 
sufficiently  novel  to  include  in  this  Report.  Here,  we  derive  equations  involving  stress  resultants  and  stress 
couples  from  the  fundamental  equilibrium  equations  (8.58)-(8-60),  which  we  repeat  here  for  convenience: 


S(0)RZi,Z  + 


i  -  S(0)ee«  +  5(o)fi©*.©]  -  °> 


(8.122) 


S(o)ezi,z  +  S(o)nei),R  +  ^  5'(o)eei,e  0- 

-j^-S(0)ezi  +  S{0)zzx]  z  +  ^  [w,flS(o)Rei  +  ^5(o)eei  +  s(o)ezi]  e 

+  ^  [R  (u,RS(0)RRi  +  ^-5(1  O)R0i  +  5(0)RZi)]  R  +  Poi  9  =  0, 


(8.123) 

(8.124) 


where  the  last  equation  is  a  slightly  modified  version  of  (8.60).  Each  of  equations  (8.122)-(8.124)  is  first 
integrated  through  the  thickness  to  eliminate  the  Z-dependence.  Introducing  radial,  circumferential,  and 
in-plane  shear  stress  resultants,  defined  by  the  following  integrals  through  the  thickness  (refer  to  the  lower 
portion  of  Figure  2): 


as  well  as 


N“-  J 

rh/1 

/  S(0)RRidZ  - 

-h/2  •> 

Ahc-h.)/2 

/  S(0)rrc  dZ  +  l 

’-h/2  J 

,h/2 

S(0)RRe  dZ, 

(hc— h»)/2 

(8.125) 

Nr6  = 

II 

§ 

<E 

* 

O 

to 

t _ 1l 

,(hc-h.)/2 

/  S(0)RecdZ  + 

J —h/2 

fh/2 

/  S(fi)RQsdZ, 

/(h„— h.)/2 

(8.126) 

Ne  = 

J 

rfc/2 

/  5'(o)e©i  dZ  — 

’—h/2 

r(hc—h,)/2 

/  5(0)e©c  dZ  +  , 

/-h/2  J 

rh/2 

/  5(0)60*  dZ, 

(he— h.)/2 

(8.127) 

out-of-plane  shear  stress  resultants  defined  by 

Qr  = 

/■h/2 

/  S(0)RZidZ  = 

/-h/2 

r(hc—h,)/2 

/  S(0)RZc  dZ  + 

/-h/2  J 

fh/2 

/  5(o)R2r8  dZ, 

\hc-h,)/2 

(8.128) 
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S(o)ezi  dZ 


f{hc 

S(o)ezc  dZ  +  I  Si0jezs  dZ, 

-h/2  J(hc-ha)/2 


equations  (8.122)-(8.124)  reduce  by  this  integration  process  (after  multiplying  thru  by  R)  to 

(RNr)  r  -  N&  +  (Vr©,©  =  0, 

(R2Nr&)  r  +  RNe,e  =  0, 


(8.129) 

(8.130) 

(8.131) 


pR  +  j R  (uJ'RNr  +  ~~ Nrq  +  Qr  +  (w^rNrq  +  -jjp-N©  4-  Qe'j  +  jogR  =  0,  (8.132) 

where  we  have  also  applied  the  leading  order  boundary  conditions  of  pressure  (8.52).  In  the  last  equation 
we  have  introduced  the  pressure  difference  p,  defined  earlier  in  equation  (8.75),  and  have  assumed  the  mass 
densities  to  be  constant  through  their  respective  thicknesses  (although  they,  as  well  as  p,  may  vary  with  R 
and  0). 

We  next  multiply  equations  (8.122)  and  (8.123)  by  Z,  to  obtain 

Z  {S(o)RZi)  z  +  ^  [(i725(0)fiRi)ifl  -  Z  S(0)eei  +  (Z  S(o)mi) ,©] 

=  (ZS(0)RZi)  z  -  S(o)RZi  +  £  [  (RZ  S{0)RRi)  R  -  Z  S( o)eei  +  {Z  S(o)R©i)e]  =  0, 

(8.133) 


Z(S(o)ezi)tZ  +  Jp  (r2  z  *%>)«©»)  ,fl  +  ^  (ZS( o)e©»)ie 


(8.134) 


{Z  S^Q-jQZi)  z  S(0)&Zi  +  jp  Z  S(o)RQi)  R  +  ^  (z  S(0)&&i)  q  —  0. 

Integrating  each  of  these  equations  through  the  thickness,  applying  the  boundary  conditions  of  pressure 


again,  and  introducing  the  following  stress  couples: 

Mr  = 

J 

/■h/2 

/  Z  S(o)  Rm  dZ  = 

’  —h/2  J 

f(hc-h,)/2 

/  Z  S(p)RRC  dZ  +  J 

'-h/2  J 

rh/2 

'  Z  5(o)rrs  dZ, 

(hc-h.)/2 

(8.135) 

Mr@  = 

fh/2 

/  Z  S(0)R©f  dZ  = 

J-h/2 

Ahc—hs)f2 

/  Z  S(0)ii©c  dZ  + 

J-h/2 

rh/2 

/  Z  S(O)R0S  dZ, 

l(hc-h.)/2 

(8.136) 

Me  = 

ph/2 

/  Z  5(o)eef  dZ  = 

' —h/2  j 

f(hc-ha)f  2 

/  Z  5(o)e©c  dZ  +  i 

f~h/2  J 

M2 

1  Z  S(o)ee*  dZ, 

(hc-h,)/ 2 

(8.137) 

we  obtain  from  (8.133)  and  (8.134)  the  following  equations  involving  the  shear  stress  resultants  and  couples: 

-Qr  +  [(rMr)r  -  M©  +  Mr©^]  =  0,  (8.138) 

and 

-Qe  +  ( R 2  Mr@)  R  +  4  M©  ©  =  0.  (8.139) 


-Qe  +  (I?2  Mr@)  r  +  j^Mete  =  0. 
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The  last  two  equations  can  be  used  to  eliminate  Qr  and  Q©  in  equation  (8.132),  to  bring  it  to  the  form 
[ii  (uj^Nr  +  -j^-NRe')  +  ( RMr)  R  -  Me  +  MRe,e J  R 


+  W'RNrb  +  ~p~Ne  +  d2  (-R2  Mrs)  r  +  »Me,e  +  (p  +  7off)B  —  0 


R 


R? 


R 


(8.140) 


,e 


Equations  (8.130),  (8.131),  and  (8.140)  are  the  fundamental  equilibrium  equations  in  terms  of  stress 
resultants  and  couples,  which  are  given  by  the  integrals  (8.125)-(8.127)  and  (8.135)-(8.137).  Performing 


these  integrals,  we  obtain  (see  Appendix  A  for  the  details): 

Nr  =  Af  +  Ae$j  R  +  Av  £q0  +  BJsrr  +  Buke&,  (8.141) 

Nrq  =  A©  e°R0  +  Be^Re,  (8.142) 

Ne  =  A/"  +  Av  (?rr  +  Acqq  +  ByliRR  +  BA:©©,  (8.143) 

Mr  —  -M  -  BeRR  -  Bve%e  -  DkRR  -  D„hee,  (8.144) 

Mrs  =  —Bee°Re  ~  Dq^rb,  (8.145) 

Me  =  —M  -  B„e°RR  -  Be%Q  -  D„kRR  -  Dkee,  (8.146) 

where  we  have  introduced  the  following  new  constant  (having  units  of  N  •  m/m): 

M  =  ±hch,(Sc  -  St),  (8.147) 


all  other  constants  having  been  defined  in  (8.116)-(8.121).  Substitution  of  these  expressions  in  equations 
(8.130),  (8.131),  and  (8.140)  yields  equations  (8.113)-(8.115),  which  were  obtained  in  the  previous  Subsection 
as  continuity  conditions. 


9  Geometrically  Nonlinear  Membrane  Shell  Laminate  Theory 


We  follow  Erbay  [11],  and  set  the  scaling  exponents  of  Section  7,  Table  1,  to  r  =  m  =  1/2,  l  =  n  =  1,  t  =  3/2, 
p  =  2,  q  =  5/2,  and  p  =  1/2,  to  obtain  a  generalization  of  the  geometrically  nonlinear  membrane  theory 
of  Hencky  [13]  and  Campbell  [14]  to  a  membrane  laminate.  With  this  choice  of  exponents,  the  constitutive 
equation  (7.23)  takes  the  form 


e5'22iSiZz  =  eS<Sr  +  £  j  (1  -  v>)  e"1'2  W,c  +  \e~l  W2  + 


+  Vi 


'U?,  +  VZ  .  (fte-^)2  +  (f.e  +  fr): 

2  +  2  p2 


)])■ 


(9.1) 


Each  variable  is  now  expanded  in  an  asymptotic  series  in  powers  of  s1/2  (recall  that  p  =  1/2  here),  and  in 
particular  we  have  for  the  partial  derivatives  of  W  with  respect  to 

Wx  =  W(0),c  +  e1/2W(iU  +  eW(2),c  +  e3/2W(3)>c  +  e2W{i)  +  0(e5/2),  .  (9.2) 
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hence,  to  order  e2,  we  find 

=  W(0U  +  2%,C%)^1/2  +  [Hu  +  2%)X%u]  £ 

+  2  [%),<%),C  +  %),(%,(]  e3/2  +  [4+2%%  +  2%,s%,J  £2  +  °(£5/2)> 

(9.3) 

and  similar  expressions  for  U^-  and  V2.  Substituting  in  (9.1),  we  obtain  to  order  e: 

£5^2£tS(o)jZ2  =  ehiS™  +  5i|(l-i/j)  -e  1  ^  +  £-1^2  +  W(o),c^(i),<) 

+  \  (2  %H  +  Hu  +  2  %U  %)X  +  UfoH  +  9<ki  ) 

+  el/2  ( W(2),<  +  W(0),c  Wfs), (  +  M7(i),c  W(2),<  +  (7(o),c^(i),c  +  ^(0),C^(1),C  ) 

+  —£  (2  W(3),<  +  Wf2),c  +  2W(i))C1¥(3)>c  +  UfV)i  +  2U(ouU(2U  +  V(\)x  +  2V{0)XV{2)X  ) 

+  «  (^+^,+3=“^  +  %*)  }  + 

(9.4) 

Since  this  expression  must  hold  in  the  limit  e-iO  (vanishing  thickness),  the  coefficients  of  negative  powers 
of  e  must  vanish,  hence  =  0,  yielding  the  solution 

%)  =  w(p,0),  (9.5) 

where  w  is  an  arbitrary  function  of  p  and  0  only.  Equation  (9.4)  then  reduces  to 

e5/2^S(0)izz  =  eXiS™  +  £  j  (1  -  **)  \  (2 %U  +  Wfax  +  %.c  +  9k< ) 

+  e1/2  (  W( 2),<  +  Wm  W(2)X  +  U(0)xU(i)x  +  V(o),cV(i),c  ) 

+  —£  (2  W(3),C  +  W(2),C  +  2W(i)iCT¥(3))C  +  %u  +  2t/(0)iCt/(2),C  +  i/(21)iC  +  2V(0)lCP(2),(  ) 


I  v-r(fj  1  1\y2  I  V(0).©  +  f/(0)  ,  W(fi),e 

+  1  l  UW’»  +  2  WW+  +  p  +  ~2^~ 


+  0(e3/2). 


The  following  three  additional  conditions  are  obtained  from  the  vanishing  of  the  coefficients  of  e°,  e1/2,  and 
£,  respectively: 

2W(i),c  +  w(i)x  +  ^(o),c  +  ^(o),c  =  °>  (9-7) 

W( 2),C  +  W(i),<  ^(2),<  +  ^(o),c^(i),<  +  ^(o),c^(i),<  =  0,  (9.8) 

2  EiS"m  4-  £i  (1  -  ^  2  W7(3),c  +  W(2),<  +  2  W(i),c  W(3))C  +  Uf1)x  +  2  f/( o),<^(2),<  +  F(2)  c  +  2V(0),cV(2),{) 

0„.  ( f>.  .  .  1..2  ,  %),e  +  ^(0)  ,  „ 


+  2  Vi  \  (7(o) +  -tu)P  + 
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which  we  shall  need  later.  We  turn  next  to  the  two  off-diagonal  constitutive  relations  (7.24)  and  (7.25)  under 
these  new  scalings: 


e2S iSRzi  =  y  [0,(  +  +  f ,PW,(  +  W,pWt<  +  e ( U,PU,<  +  W C )] . 


22 iSezi  =  y  | Vx+e1/2^ 


+  Me+£ 

P 


VAU.e-V)  +  Vx<y,e  +  V) 


]} 


(9.10) 

(9.11) 


The  leading  order  terms  in  these  equations  imply  (7(o),c  =  0  and  ^0),<  =  0,  hence  the  leading  order  in-plane 
displacement  components  are  both  independent  of  C,  i-e-i 

U{0)  =  u{p,@),  V(o)  =  v(p,G). 

The  conditions  C/(0),<  =  0  and  7(0),c  =  0  reduce  equations  (9.7)-(9.9)  to 

(2  +  %u)  =  0. 

%U  (i  +  %u)  =  o, 


2  T,iS?m  +  Si 


(1  -  Vi)  (2W(3),(  +  W?2U  +  2  W(i)>cW(3),<  +  +  ^(i),c) 


+  2vi(u,p  +  \wl  +  ^y^  +  ^j  =  0. 


(9.12) 

(9.13) 

(9.14) 

(9.15) 


Noting  that  W(0),?  =  0  implies  a  scaled  Jacobian  determinant  (7.22)  of  the  form 

J  =  i  +  %u  +  °(£l/ 2)> 

it  follows  that  in  the  limit  e  ->■  0,  since  we  must  have  J  >  0  in  that  limit,  we  must  also  require  WliC  >  -1. 
This  precludes  the  solution  W(i)iC  =  -2  of  (9.13),  implying  that  W{1)>(  =  0.  This  leads  to  W^)x  =  0,  from 
(9.14),  hence  (9.15)  yields  the  following  expression  for  W(3)iC: 

-  -5 (%,  +  %-c)  -  ^  k  ♦"?=+£)  •  t9M) 

Next,  the  off-diagonal  in-plane  constitutive  relation  (7.26)  has  the  following  form  under  these  scalings: 

U,e-V\  .  JJr,e  .  W,pW,b 

(9.17) 


eXiSRQi  =  y  |«  ^  +  G'&p  -) 


+  e  g  +  e- 

P  P 


+  er 


UAU.e-V)  +  VAV,e  +  V) 


])- 


The  leading  order  term  here  yields  the  same  constitutive  relation  found  in  (8.18),  except  that  the  strain- 
displacement  relation  has  changed: 


Gi 


Sj5(o)n©i  =  y  l  v>p  + 


u,e~v  ,  r ,pu\e  w,Pw,e  \  _  r 


+  + 

P  P  P 


=  Gi  ej?e, 


(9.18) 
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where 


?  _  1  [  „  u>e-v  T<pw,e  w,Pw,e\ 

2  V  p  p  p  ) 

The  final  two  constitutive  relations  (7,27)  and  (7,28)  take  the  following  forms: 


eXiSRRi  =  eT^ST  +  ft  (1  -  **>  i  ^  +  sTJ¥iP  +  -  ^  +  V; 


+  v%  <  £ 


P  /  2p2 


(?7,e-y)2  +  (^e  +  ^) 

2p2 


+  +  t-iwj  +  j(5.<+^)})’ 


eSiSeei  =  eEiST*  +  ft  (1  -  *«)  «|  e 


£e+gV  c^,2e  t  .2  [ (£e - ^)2  +  (ge  + 1/)2 


+  **  |  +  l  +  e2  (  u%  +  V2P  )] 

+  *~1/2^C  +  +  \  (V< +  K<)  }  V 


In  the  last  two  equations  the  terms  up  to  0(e)  involving  derivatives  with  respect  to  £  reduce  to  simply 
W( 3}^j  which  can  be  replaced  by  (9.16).  The  coefficients  of  e  then  yield,  after  some  algebra,  the  same 
in-plane  constitutive  relations  (8,27)  and  (8,28),  except  that  the  strains  are  now  given  by  the  ^-independent 


expressions 


tRR  =  ,  €©e  = 


vjb±u  w% 


The  most  important  results  here  are  that  the  leading  order  in-plane  displacement  components  have  the  forms 
U(q)r  =  u(R,  0)  and  t/(o)©  =  v(Ry  ©)  (in  terms  of  the  physical  variables),  i.e,,  they  are  independent  of  Z . 
This  is  in  marked  contrast  to  the  KirchhofF-Love  expressions  (8.31)  found  previously,  which  are  linear  in  Z. 
As  a  result  of  these  simple  expressions  for  the  in-plane  displacements,  the  curvature  terms  kRR,  and 
k@e  do  not  appear  in  equations  (8.35)-(8.37)  defining  the  in-plane  strain  components.  Summarizing  to  this 
point,  we  have  derived  the  following  expressions  from  the  scaled  constitutive  relations: 

U(o)z  =  w(RyB),  U(0)R  =  u(J?,0),  l7(o)e  =  v(i?,0),  (9,23) 

S(o)RQi  =  GitRs ,  (9,24) 

E* 

S(0)RRi  =  Si  +  Y~^/2  (  £RR  ^i6©©  )  >  (9-25) 


5(O)©0i  =  Si  +  (  £©©  +  )  : 


where  S»  is  given  by  equation  (8.29),  and 


If  ,  «,©  -  v  ,  ,  w„R«>,e 

eiie  =  5  K*  +  — r—  +  — b“  +  — r— 


-  r° 

=  efi©i 
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(9.28) 


tRR  =  U.fl  +  r  ,RWtR  +  ~W2r  =  e°RR, 

£ee  =  R  +  2&  =  £ee-  (9,29) 

Under  these  new  scalings,  the  radial  and  circumferential  equilibrium  equations  (8.122)  and  (8.123)  are 
ichanged,  but  the  shear  stresses  disappear  from  the  terms  of  the  axial  equilibrium  equation  (8.124)  involving 
irivA.t.ivps  with  rfisnect.  to  R  and  A.  The  leading  order  eauilibrium  eauations  are  then  riven  bv 


S(o)RZi,z  +  jj;  [(-R5(o)hr,)  R  -  5(o)©ei  +  <S(o)fle«,e]  —  0, 

(9.30) 

S(o)ezi,z  +  jp  ^  S(o)nei),R  +  j j  S(0)ee»,e  =  0. 

(9.31) 

(v,rS(o)rzi  +  —j^-S(0)ezi  +  S(o)zzi}  ^  +  —  (v,nS(o)iiei  +  -jp-S(o)eei)  Q 

1  r  /  W  9  V 

+  S  l  y°’RS(o )RRi +  )Rei ) 

1  +  Poi  9  =  0,  (9.32) 

J  ,R 

where  we  used  the  result  +  0(e2)  in  obtaining  the  scaled  version  of  the  last  equation.  The 

boundary  conditions  of  pressure  are  again  given  by  (8.53),  so  the  equilibrium  equations  can  be  integrated 
through  the  thickness,  applying  the  boundary  conditions  of  pressure,  to  obtain  (after  multiplying  thru  by 
R)  equations  in  terms  of  the  in-plane  stress  resultants  (8.125)-(8.127): 

(RNr)  r  -  Ne  +  NRe,e  =  0, 


{R2NRe)R  +  RNe,e  =  0, 


(9.33) 

(9.34) 

[■ R  (u,rNr  +  R  +  (u,rNRq  +  -^jplVe)  0  +  (P  +  1o9)R  =  0.  (9.35) 

Since  the  curvatures  have  vanished  from  the  formulation,  the  stress  resultants  (8.141)-(8.143)  in  this  case 
reduce  to 

Nr  =  N  +  Ae®RR  +  Aue%&,  (9.36) 

NRe  =  Ae  e°Re,  (9.37) 

Are  =  A/"  +  Av^rr  +  A  £00 .  (9.38) 

10  Geometrically  Linear  Shell  Laminate  Theory 

Assuming  the  scaling  exponents  r  =  l,  m  =  2,  £  =  n  =  3,  t  =  p  =  4,  q  =  5,  and  n  =  1,  we  begin  again 
with  the  constitutive  relation  (5.5)  for  Szz,  from  which  we  obtain 

esE iSzzi  =  e^iS™  +  £i{  (1  -  Vi)  eW,(  4-  e2^  +  e4  ^  (t%  +  V2(  ) 

Wl 


",  [«3  (  o,,  +  r,fw,„  +  +  f*  (f  i V%  +  ^ 


+  +  )]}. 


( 


+  n  .  (U,e  -  U)2  +  (Ue  +  V)2 

2p2 


(10.1) 
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Each  of  the  scaled  stress  tensor  components,  as  well  as  each  of  the  displacement  components  and  its  partial 
derivatives,  is  now  written  as  an  asymptotic  expansion  in  S  =  e,  since  /i  =  1.  Substituting  the  appropriate 
expansions  for  the  variables  appearing  in  equation  (10.1),  we  obtain  to  third  order  on  the  right-hand  side: 


e5ViSmzi  =  £3Si5fm  +  £i{(  1  -  **) 


eW(o)6  +  **  (2 + 


+ 


(i  -  vi)  (w{2U  +  w(ou  wau)  +  Vi  (um>p  +  fpw(0)tP  +  Ym±Bsi 


+ 


0(e4) }  H0.2) 


The  leading  order  term  of  this  relation  yields 

W{ou  =  0  =>  W(  0)  =  w  (p,0),  (10.3) 

’where  w  is  an  arbitrary  function  of  p  and  0  only,  and  then  the  next  order  term  reduces  to  the  result 
=  0.  With  these  results,  the  third  order  term  then  gives  the  following  important  relation: 

Vi 


rxr  onm  tt  \  V  in  I 

-  -£i(i-vjSi  ~  [Uw’p + T'pw'» +  p 


(10.4) 


Next,  consider  the  last  two  off-diagonal  constitutive  relations  of  (5.6)  for  the  out-of-plane  stress  com¬ 
ponents.  The  scaling  of  these  components  yields  the  scaled  relations  (to  third  order  on  the  right-hand 
sides): 


e^iSnzi  =  y  [e2  (%),<  + W{0U  +  T,PW{ m)  +  0(e3)] , 


(10.5) 


e^iSezi  =  $ 


s2  I  V, 


W, 


(o),C  ' 


(o),e 


+  0(e3) 


(10.6) 


respectively.  From  the  leading  order  terms  of  these  expressions,  recalling  that  W(o),£  —  0,  we  obtain  the 
following  two  equations: 

W  Pi 

U(0U  +  ™,p  -  0,  V(0).<  4-  -y  =  0, 


which  may  be  integrated  to  obtain 


w,e 


U( o)  =  u(p,  0)  -  C  w,p  ,  and  V(0)  =  v(p,  0)  -  C 

r 

where  u  and  v  are  arbitrary  functions  of  p  and  0  only. 

The  scaled  version  of  the  first  off-diagonal  in-plane  constitutive  relation  of  (5.6)  reduces  to 


e3£jSnei  =  Gi 


V +  Emjt^rn  +  EAe 

i  P  P 


+  0(£4) 


Thus,  to  leading  order  we  obtain  the  simple  relation 

£ iS(0)RBi  =  GitRe , 

after  introducing  an  in-plane  strain  component  ?b©  defined  by 

-  _  1  (c,  ,  ^(0},e  -  f  ,pW,e 

ene  =  ^  I  Uo),p  +  - - - +  — ~ 


1 


+  +  E>£  +  2C 

P  P 


( _  w>Pe\ 

\f ?  P  ) 


(10.7) 


(10.8) 


(10.9) 
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The  final  two  constitutive  relations  (5.3)  and  (5.4)  take  the  following  forms  (to  third  order  in  e),  after 
scaling  of  the  various  components: 

esXiSRRi  =  e3XiSr  +  £<{(1  -Vi)e3  (t/(0),„  +  T,pWWiP) 

+  *  [e3  (  +  £  (w(ou  +  eWwx  +  e2%u) 

+  \e2  +  2eW(0),cW(1),<)  ]  +  0(e4)  },  (10.10) 

and 


+ 

+ 


e3Ej5eei  =  e3ZiS?m  +  6  {  (1  -  »i)  e3  ^  ^(0)’9  j 

Vi[e3  ( l/(o),  P  +  +  e  (W(0),c  +  eW{  i),<  +  £2^(2),<) 

5«*  (w?o),C  +  2eWmW(1)x)  ]  +  0(e4)  }. 


(10.11) 


The  first  and  second-order  terms  involving  in  these  expressions  vanish  since  W(o),<  =  W(i)^  =  0,  and 
(10.4)  can  be  used  to  replace  the  third  order  term  involving  To  simplify  the  reduction  of  the  leading 

order  expressions  that  follow,  we  introduce  in-plane  strain  components  crr  and  €©e  defined  by 

€rr  =  JJ(  o),p  4-  TiPwp  =  2)P  +  TiPwyP  —  Cw,pp  =  —  C^rr,  (10.12) 


%  fkejJk  =  £e±“  _  c  (fit  +  %)  =  a«  9  -  cfee,  (10.13) 

P  P  \  P  P  ) 


where  we  have  introduced  ^-independent  terms 

f$iR  —  u,p  +  r  ,pW,p,  kRR  =  w<pp, 


w)  _  v,e  +  u 

e©e  =  — —  - 


P  P 


These,  together  with  (10.4),  allow  us  to  write  the  leading  order  relations  as 

XiS(0)RRi  =  XiS?m  +  Si  [(1  -  Vi)eR R  +  vi?ee  -  (?RR  +  «©e) 

s  (l^)  *lSr  +  £t  +  ^  ■ 

Si5(o)ee»  =  (jr£)  m  +  Si  [(1  -  i/<)?ee  +  VitRR  -  ( eRR  +  ?ee) 

s  (t^t)  SlS?” + *  (t^t)  + ■'«**>  ■ 

Replacing  Si  in  the  last  two  expressions  yields 

^  ^  E' 

£ iS(o)RRi  =  £»■£*  +  _*  g  (?rr  +  Vitee ) , 

A  Vi 


(10.14) 

(10.15) 


(10.16) 
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Et 


EiS( o)©04  —  Sj  Si  +  "2  (  £0©  +  Vi^RR  )  . 

* 


where  we  have  introduced  the  new  constants 


*  s  (££)  *-• 


(10.17) 


(10.18) 


We  summarize  to  this  point  by  rewriting  the  important  results  in  terms  of  leading  order  variables  that 
are  functions  of  the  physical  coordinates  R,  0,  and  Z,  viz.. 


u  =  e3au,  v  =  e3av,  w  =  e2aw,  T  =  ear, 


(10.19) 


S(o)R@i  =  e3SfS( 0)ii©i  =  Gi  £r©- 


(10.20) 


S(o)RRi  =  £3Si5(0)j?fli  =  Sj  +  +  ^£ee ) > 


(10.21) 


where 


■5(0)001  =  £  SjS(0)©©i  =  Si  +  -  —  g  ( £0©  +  vitRR ) 


(10.22) 


(10.23) 


and  the  last  equality  of  (10.23)  holds  for  the  special  case  of  thermal  stress^  following  from  (5.2).  The  strain 
components  are  given  by 


3^v.  l 

£r&  =  £  £r®  —  2 


..  ,  ue-v  TRwiB  t  nr7  /wye  wtRe\ 

v-R+—ir+—R-*iZ\-w-—) 


=  e%e  -  ZkRe,  (10.24) 


CflR  =£3tRR  =  UtR  +  VtRWtR  —  Zw%RR  —  €%r  -  ZkRR, 


(10.25) 


_  3-~-  V,  ©  +  u 

€@©  =  e  e@©  =  — — — 


where  the  S-independent  strains  and  “curvatures”  are 

,o  _  1  ( ..  ,  u,e~v  T,RWte\  ,  _  w,e  ,  to, a© 
e*e  =  ^  \V’R  +  b~  +  ~R~ )  ’  k™  =  “1^  +  “r". 


(10.26) 


(10.27) 


eRR  =  U,R  +  r,RW,R  ,  kRR  =  W'RR, 


(10.28) 


„o  _  v,&  + u  ,  -  W,R  ,  w,ee 

e&e  =  ^~’  kee  =  ^  +  ~W 


(10.29) 


Applying  the  scaling  exponent  values  of  this  Section  to  the  scaled  equilibrium  equations  (7.38)-(7.40), 
we  obtain  to  leading  order: 


■5(0 )RRi,p  +  -S(Q)RQi,e  +  S(0)RZi,(  +  ~  (§{0)RRi  ~  S(0)©©2)  =  0, 


(10.30) 
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(10.31) 


—  1  2  ^ 

S(0)RQi,p  +  “-5(0)601,6  +  S(0)62i,C  +  “5(o)fl©i  =  0, 


(r ,pS(0)RRi  +  S(0)RZi)  +  -  (r,pS(0)flei  +  5(0)6 Zi) 

+  (r,p5(0)flZi  +  5(0)ZZt)  +  -  fr,pS(0)RRi  +  5(o)flZi)  +  9  =  0.(10.32) 
In  terms  of  the  physical  variables,  these  take  the  forms 

S(0)RZi,Z  +  ^  [(-R5(0)fifii)  fl  -  5(o)©©t  +  5(o)fl©j,©|  =  0,  (10.33) 


S(0)6Zi,Z  +  Jp  ^(0)Uei),H  +  ^  5(0)©©i,©  =  0. 


(10.34) 


[r,ftS(0)flZi  +  5(0) zzi]  z  +  ^  [r,j?5(o)fl©i  +  5(0)©zi]  e 


+  R  (^'fl5(°)flfli  +  5(°)flZt )]  R  +  Poi9  =  0,  (10.35) 

The  leading  order  boundary  conditions  of  pressure  are  the  same  as  for  the  previous  theories,  viz., 

^(o )rz  =  0’  £(0 )ez  =  S%)zz  ~  ~P±-  (10.36) 

Integrating  equations  (10.33)-(10.35)  through  the  thickness,  we  obtain  the  equilibrium  equations  by  the  same 
method  used  in  Section  8  in  terms  of  the  stress  resultants  and  couples  introduced  in  equations  (8.125)-(8.127), 
and  (8.135)-(8.137),  viz.: 


(RNR)tR  —  Nq  +  Nrs,b  =  0, 
+  RNe,e  =  0, 


(10.37) 


(10.38) 


and 


\rT,r  Nr  +  ( RMr)  R  -  Me  +  AO*©,©]  r 

+ 


r ,rNrs  +  jp  (R2  Mrs)  r  +  -^M©,© 


J  ,© 


R?  v . ™'’R  ■  R 

The  stress  resultants  and  couples  are  again  given  by  equations  (8.141)-(8.143),  i.e., 
Nr  =  Af  +  AeRR  +  Av^qq  +  B  kRR  +  Bpkee, 

Nrb  =  Aq  eRe  +  BekRe, 

Ne  =  N  +  A„  eRR  +  A  eg©  +  Bv  Urr  +  B  kee  , 

Mr  =  -M  -  Be°RR  -  Bp  £00  -  DkRR  -  Dvkee, 


+  (p  +  7o  9)R  =  0.  (10.39) 


(10.40) 

(10.41) 

(10.42) 

(10.43) 
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(10.44) 


Mrs  —  -Be  e^e  -  DekRe, 


Me  —  —  M  —  Bve°RR  -  Be%e  —  DvkRR  -  Dkee,  (10.45) 

but  note  that  the  strains  and  curvatures  are  now  replaced  by  (10.27)-(10.29)  of  the  present  Section. 

11  Geometrically  Linear  Membrane  Shell  Laminate  Theory 

In  this  Section,  we  develop  the  equilibrium  equations  for  a  geometrically  linear  membrane  laminate, 
subject  to  pressure  and  gravity  loading.  We  begin  by  choosing  the  scaling  exponents  of  Section  7  to  be 
r  —  m~l,  i  =  n  —  3/2,  t  =  5/2,  p  =  3,  q  =  7/2.  and  p  —  1/2,  as  indicated  in  Table  1  at  the  end  of  that 
Section.  With  these  exponents,  equation  (7.23)  takes  the  form 


£7/2ZiSZzi  =  e3/2£  iS?m  +  £ 


Vi) 


w,c  +  \K  +  ‘\{vi  +  v}) 


+  -<  k/!  (u,f  +  +  e!  [h^„  +  H  'j 


(11.1) 


U%  +  K  ,  (Ee^Vl+JXe  +  Uf 

2  2p2 


The  terms  involving  partial  derivatives  with  respect  to  C  must  be  expanded  up  to  terms  of  order  £3/2,  as  in 
equations  (9.2)  and  (9.3).  Making  the  appropriate  substitutions  in  (11,1),  and  then  setting  the  coefficients 
of  e°,  £1,/'2 ,  e,  and  £3;2  in  equation  (11.1)  to  zero,  we  obtain  the  following  four  equations: 


"W  +  2^fo),  c  =  °- 


W(1),C  +  W(0),<;  W(i),c  =  o, 


^(2),C  +  ^(21),C+2^(0),C^(2),C)  +  (qkc  +  ^ou)  = 


(11.2) 

(11.3) 

(11.4) 


W  +  £i|  (1  -  Vi)  JW(3),<  +  (  W(0)xW(3)£  +  W(l),<;  ll/(2),c ) 

+  (p(o)x  U(1)X  +  V(0)X  V(1)>C  )  ]  +  Vi  ^(0),p  +  -(0),e^(0))  |=0.  (11.5) 

Equation  (11.2)  implies  that  W{0)iC  =  0,  hence  W{0)  =  w(p,  Q,  and  (11.3)  then  yields  t?(1)iC  =  0,  reducing 
the  remaining  two  equations  to 


W(2),C  +  (^(0),c  +  ?o),<  )  =  0, 


(11-6) 


SiSr+^|(i-Vi)[w(3)lC  +  (&(0),ct/(1))<  +  +  |=0.  (11.7) 

We  turn  next  to  the  two  off-diagonal  constitutive  relations  (7.24)  and  (7.25)  under  these  new  scalings: 

e^SRZi  =  ^  [el>2  Ux  +  eWp  +  eftPWx  +  eWpWx  +  e2  ( U,PUX  +  VPVX  )]  ,  (11.8) 
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e3Si5ez 


■  -!{ 


P  P 


Ux&e-V)  +  VAV>e  +  U) 


(11.9) 


The  leading  order  terms  in  these  equations  imply  E/(o),c  =  ®  an<^  ^(0),c  =  hence  the  leading  order  in-plane 
displacement  components  are  both  independent  of  £,  i.e.,  U^q)  =  u(p>®)  and  ^(o)  =  w(p>©)*  When  these 
two  results  are  substituted  in  (11.6)  and  (11.7)  we  find  that  W(2)x  =  0,  and  is  given  by 


^  s Pi  (  tt 

=  '60^  "  v(0)'' 


^(o),e  +  U(o) 


)■ 


(11.10) 


which  will  be  used  later.  First,  we  note  that  the  off-diagonal  in-plane  constitutive  relation  (7.26)  has  the 
following  form  under  these  scalings: 


e3/2ZiSRe, 


_  Gif 
2  1 


U.e  -  V\  .  2  rtfiW,e  . 


e3/2  o  +£: 


+  ea 


+  e 

P 

uMe^]  +  VJ^e±U) 

P 


11 


(11.11) 


The  leading  order  term  here  yields  the  same  constitutive  relation  found  in  (8.18),  except  that  the  strain- 
displacement  relation  has  changed: 


£<5(O)rt0i  =  "y  (v,p  +  p  =  Gi 


£rq, 


€rs  =  2  ^  + 


17,  e  ~  V 


(11.12) 


The  final  two  constitutive  relations  (7.27)  and  (7.28)  take  the  following  forms: 
es/*XiSRRi  =  e^XiSr  +  £  ( (1  -  *i)  {  ^  U,p  +  e2f ,PW,P  +  \  [e2W2p  +  e3  ( U2p  +  V2  )] } 


/  -.3/2  (V,e  +  U\ 

'(U,e-V)*  +  (V,e  +  U)2] 

{  \  >  ) 

1  +£  2?  +£ 

1 

to 

(11.13) 


+  WA  +  ^  W%  +  £ 


^ +?<)})• 


and 


„  W%  .  , 

+  £  V+£ 


e3/2S iSeei  =  e3/2£iSfm  +  6  ^(1  -  *)  je3/2  (^7^) 

+  *J e3/2 U,p  +  e2f ,PW,P  +  \  [e2W2p  +  e3  ( U2p  +  V2 )] 

+  W,(  +  \w2+el-(u2<:  +  V2)^y 


(ge  -  V?  +  (Y,e  ±  U? 
2  P2 


(11.14) 


The  terms  involving  derivatives  with  respect  to  £  in  the  last  two  equations  reduce  to  simply  W(3),C’  which  can 
be  replaced  by  (11.10).  The  coefficients  of  e3/2  then  yield,  after  some  algebra,  the  same  in-plane  constitutive 
relations  (8.27)  and  (8.28),  except  that  the  strains  are  now  given  by  the  simpler  expressions 


£rr  =  UiP,  eee  = 


V,e  +  U 


(11.15) 
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As  with  the  geometrically  nonlinear  membrane  theory,  we  again  find  that  the  leading  order  in-plane 
displacement  components  ?7(o)k  =  u(R ,  0)  and  I/(0}©  =  v(R,  0)  are  independent  of  Z.  The  curvature  terms 
kRR ,  kue,  and  fc@©  are  thus  again  absent  in  equations  (8.35)-(8.37)  defining  the  in-plane  strain  components. 
In  addition,  the  nonlinear  terms  involving  derivatives  of  w,  as  well  as  any  reference  to  the  surface-defining 
function  T(R),  have  disappeared  from  the  strain-displacement  relations.  In  terms  of  physical  variables,  the 
governing  equations  can  be  summarized  thus  far  by: 


U(o)z  =  w(R,Q),  U(0)R  =  u(R,  0),  Z7(o)e  =  «(#,©)> 

(11.16) 

S(o)Rei  =  Gi  eRQ, 

(11.17) 

E’ 

S(o)RRi  =  Si  +  j  ( eRR  +  I'ieee ) , 

(11.18) 

E- 

S(o)e&i  =  Si  +  2  ( €ee  +  VitRR ) , 

(11.19) 

where 

If  t  u^e-v\  _  0 

efle  =  2  1  v, R  +  R  1  -  eRe, 

(11.20) 

CRR  =  u,r  =  eRR, 

(11.21) 

v,e  +  u  _  o 

€&B  ~~  JI  —  eBB‘ 

(11.22) 

Under  these  new  scalings,  the  radial  and  circumferential  equilibrium  equations  (8.122)  and  (8.123)  are 
unchanged,  but  the  out-of-plane  shear  stresses  disappear  from  the  axial  equilibrium  equation  (8.124).  The 
leading  order  equilibrium  equations  are  then  given  by 

S(o)RZi,z  +  ^  [(•R5(o)flflj)  Jj  -  S(0)e©i  +  ’%)H©*\e]  =  0, 

(11.23) 

S(o)ezi,z  +  -jp  (R2  S{o)Rei),R  +  -^S^ee^e  =  0, 

(11.24) 

U,RS(0)RRi  +  )R&i 

These  may  now  be  integrated  through  the  thickness,  applying  the  boundary  conditions  of  pressure,  to  obtain 
(after  multiplying  thru  by  R)  equations  in  terms  of  the  in-plane  stress  resultants  (8,125)-(8.127): 


(RNr)  R  -  Nq  +  NRete  =  0, 

(11.26) 

(R2NRq)  r  +  RNe,e  =  0, 

(11.27) 

j#  (u,rNr  + 

R  +  (u,RNRe  +  Q  +  (p  +  1fog)R  =  0. 

(11.28) 

)]  R  +  Poi9  =  0.  (11.25) 


S(o)zzi,z  + 


v,nS(o)Rei  +  —j^S(o)e&i  j  q  + 


[*( 
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Since  the  curvatures  have  vanished  from  the  formulation,  the  stress  resultants  (8.141)-(8.143)  in  this  case 
reduce  to 


Nr  =  N  +  Ae°RR  +  A^e©©, 

(11.29) 

NrQ  =  A©  6$*©, 

(11.30) 

iV©  =  M  4-  Ave°RR  +  Ae©©. 

(11.31) 

12  Conclusions 

The  method  of  asymptotic  expansions  has  been  applied  in  a  purely  formal  way  to  the  geometrically  nonlinear, 
three-dimensional  equilibrium  equations  of  a  coated  membrane  laminate,  each  component  of  which  is  assumed 
to  be  a  uniform,  homogeneous,  isotropic  elastic  material.  Four  different  choices  of  the  scaling  exponents 
appearing  in  equations  (7.2)-(7.5)  have  led  to  four  distinct  theories,  each  of  which  reduces  to  a  well-known 
theory  when  specialized  to  a  single  material.  The  method  is  systematic  and  self-consistent,  as  the  only 
freedom  available  in  constructing  a  theory  is  the  choice  of  scaling  exponents.  In  short,  having  made  a  choice, 
one  is  led  without  benefit  of  any  further  assumptions  to  a  theory  dictated  by  that  choice.  The  main  weakness 
in  the  method  is  the  essentially  ad  hoc  nature  of  the  choice  of  exponents  used  in  the  scaling  of  variables  and 
loads.  Attempts  have  been  made  in  the  recent  literature  (see,  for  example  [17,  18])  to  remedy  this  remaining 
unsatisfactory  feature. 

Although  no  attempt  has  been  made  here  to  put  the  method  on  a  rigorous  mathematical  foundation, 
articles  doing  so  are  available  in  the  literature,  especially  from  the  French  school  led  by  Ciarlet  [19,  20, 
and  references  therein]  and  his  co-workers.  Our  goal  was  to  survey,  and  systematically  derive  by  a  single 
method,  generalizations  of  several  well-known  theories  of  a  single  material  to  a  coated  membrane  laminate. 
The  mechanical  behavior  of  such  a  laminate  is  of  considerable  practical  interest  to  those  in  the  aerospace 
community  committed  to  the  manufacture  and  deployment  of  large  aperture,  optical  quality  reflectors  in 
space.  The  availability  in  a  single  publication  of  several  models  from  which  to  choose  for  the  analysis  of  such 
a  laminate  will  hopefully  be  useful  to  other  workers  in  the  community.  In  that  regard,  a  companion  Volume 
II  of  this  Report  providing  details  of  the  analytical  solutions  of  several  boundary  value  problems  associated 
with  these  models  is  currently  being  prepared  for  publication. 


Acknowledgements 

I  thank  Dr.  Dan  Segalman,  former  Program  Manager,  Structural  Mechanics,  Air  Force  Office  of  Scientific 
Research,  and  Dr.  Dean  Mook,  current  Program  Manager,  Structural  Mechanics,  Air  Force  Office  of  Scientific 
Research,  for  their  support  of  this  work. 

A  Generalization  to  a  Multilayer  Coating 

High-reflectance  optical  coatings  typically  consist  of  a  dielectric  “stack”  of  two  or  more  coating  materials 
alternating  in  position  through  the  thickness.  In  this  Appendix  we  generalize  our  work  to  include  such  mul¬ 
tilayer  coatings.  The  strain-displacement  relations  and  equilibrium  equations  characterizing  a  given  theory 
are,  except  for  the  gravitational  term  of  the  axial  equilibrium  equation,  unchanged  by  this  generalization. 
The  only  calculations  that  require  modification  are  the  through-the-thickness  integrals  of  the  constitutive 
relations,  which  define  the  stress  resultants  and  couples.  These  were  first  introduced  in  Subsection  8.6.  We 
give  details  of  the  calculations  for  Nr  and  Mr ,  noting  that  those  for  IV©,  Nrq,  M©,  and  Mrq  involve 
nothing  new,  and  will  be  seen  to  be  obvious  generalizations  of  the  earlier  single-coating  (two-layer)  results. 
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(hi ?  Vi}  Si) 


(huEl}PlfS1) 


Figure  5:  Geometry  of  multilayer  stack. 


The  geometry  of  the  reference  placement  of  our  membrane  substrate  layer  N,  with  N  —  1  coatings,  is 
illustrated  in  Figure  5  (compare  to  Figure  2),  To  do  the  through-the-thickness  integrals,  we  introduce  a 
change  of  variable  to 

£  —  Z  +  (A,l) 

As  in  Figure  2,  the  middle  plane  of  the  stack  is  the  plane  Z  —  0  (corresponding  to  £  =  hf  2)  which  has,  for 
illustrative  purposes  only,  been  placed  in  the  coating  numbered  N  —  2  in  Figure  5,  From  this  Figure,  it  can 
be  seen  that  the  thickness  hi  of  layer  i  is  given  in  terms  of  the  new  coordinate  £  by 


hi  ~  £j  -  £i— 1? 


and  the  total  thickness  is 


&  =  fa  =  -  &-i)  =  Civ- 


Thus,  for  the  through-the-thickness  integral  of  an  arbitrary  function  f(Z)  we  have: 


rh/ 2  rh  N 

/  mdz  =  /  ,/m  =  E  /  /(0«, 

Jz=-h/  2  ^€=0  i=1 


(A.4) 


where  /(£)  =  /(£  -  h/2). 

We  begin  by  making  this  change  of  variables  in  each  of  the  definitions  (8.125)-(8.127)  and  (8.135)-(8.137) 
of  the  stress  resultants  and  stress  couples,  respectively: 


/•V  2  *  rt< 

Nr  =  S(0)RRi(Z)dZ  =  E  /  S(o)RRi(Qd£, 

J—h/2  J^i-i 


(A.5) 


fh/2  N  fii 

NRe  =  /  Smei(Z)dZ  =  Y,  /  Smei(t)d(;, 

J—h/2  J €=Ci-i 


(A.6) 


/•ft/2  *  rC. 

Nq=  S(0)eei(Z)  dZ  =  E  /  ^(o  )eei(£)d£> 

J—h/2  </«=€<-, 


(A.7) 


Mr  =  /  Z S(Q,)RRi{Z)dZ  =  Y,  ($  —  h/2)  5(0)RRi 

7-h/2  j=1  ^C=€i-, 


(Odf, 


(A.8) 


/•ft/2  *  /•«- 

Mrq=  ZS(0)Rei(Z)dZ  =  Y  (Z-h/2)S{0)Rei(OdZ,  (A.9) 

J—h/2  •/«=« i-l 


/•ft/2  "  rb 

Me  =  /  ZS{0)eei(Z)dZ  =  V  /  (£-h/2)  5(0)eei 

•/-ft/2  ^7  •/?=«,-! 


(Ode- 


(A.  10) 


Recall  now  the  separation  of  each  of  the  in-plane  constitutive  relations  (8.32)-(8.34)  into  Z-dependent  and 
Z-independent  parts: 


S(0)iiez  =  CTfl©x  ~  ZrjRQii  S(0)RRi  = 

0Ri  —  ZrjRi,  5(o)ee,  = 

(A.  11) 

where  the  Z-independent  functions  axe  defined  by 

0\R0t  =  Gt^e, 

=  GikRe, 

(A.  12) 

cr/ti  =  Si  +  Qi  +  ^i€0©) , 

VRi  =  Q%  {Krr  +  Vikee) , 

(A.  13) 

<7©t  =  Si  +  Qi  (c©e  +  vie*RR) ) 

T]e %  =  Qi  (kee  +  i ^rr). 

(A.  14) 

After  the  change  of  variables  we  have 

<S(o  )RRi  =  crRi  —  (£  —  h/2)riRi, 

S(o)OGi  —  &Oi  (£  h/2)^©i, 

(A.  15) 
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and 


=  ami  ~  (£  -  h/2) rjRQi. 
Substitution  of  the  first  of  equations  (A. 15)  in  equation  (A.5)  yields 

N  r£« 


(A.  16) 


Nn  =  £[  [vm  -  (£  -  h/2)r]Ri]  d£, 

<=l  '£=£<-i 

jv  r  h  ii 

—  E  _  f*-1)  +  2^m  ~  ^i_1)  “  *?*»2  (€i  ~  &i)  ,  (A.17) 

"  r  i  1 

=  E  +  2T,Rihi^h  ~  (&  +  &-l)]j» 

where  we  made  use  of  (A,2)  to  obtain  the  final  equality.  Now,  from  (A.2)  it  is  easy  to  show  that 

i 

6  =  E  hk,  (A.  18) 

from  which 


k= 1 


i— 1 


&  +  &-i  =  E  +  E 

fc=l  fc=l 


hence 


*— l 


N 


“  (&  +  l»- 1)  =  h  -  ^2  hk  -  ^2  hk  =  ^2  hk  -  ^2  hk- 


(A.  19) 


(A.20) 


fc— 1  fc=l  fc=t  fc=l 


where  the  last  result  follows  by  substituting  the  first  equality  of  (A.3)  for  h.  Substitution  of  this  result  in 
(A. IT)  yields 


N 


N 


N 


Nr  —  &Rihi  +  53  mihi  |  53  hk  —  53  hk  . 


(A.21) 


i-l 


i- 1 


&=1 


In  this  expression  and  others  that  follow  we  make  use  of  the  following  identity,  which  can  be  proven  using 
the  principle  of  mathematical  induction: 


N 


N 


N- 1  N 


E aihi  ( E hk  ~  XI hk )  =  XI  XI  hihk  ~  °*)  > 

<  =  1  \k  =  i  k  =  1  )  i  =  l  *  =  *+! 

for  any  indexed  function  a*.  This  allows  us  to  rewrite  (A.21)  as 

N  «  N—l  N 


(A.22) 


j  if  — jl 

nR  =  53  aRi^i  +  53  hihk  {rim  -  7}Rk ) . 


(A.23) 


i=l 


t  =  1  jfe  =  i+l 


Substituting  now  for  am  and  r]Ri  from  equation  (A.  13),  we  find 

iV 

Aft  =  53  hi  [ft  +  Qi  {e%R  +  W©©)]  4- 

i= 1 

^  i¥-l  AT 

2  53  53  W*  $RR  +  ^©e)  —  Qft  {hRR  +  ^fefc©©)] ,  (A.24) 

i  —  l  £=rffl 
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which  can  be  written  in  the  same  form  as  equation  (8.141),  i.e., 


Nr  —  Af  4  A  4  Av  6@0  4*  B  fcRR  4*  Bv  &©©  , 

(A.25) 

where  the  multilayer  constants  are  given  by 

N 

v  = 

X  ^ 

i=l 

(A.26) 

N 

N 

A  =  X 

Av  =  X  hiQiVi, , 

(A.27) 

i- 1 

1=1 

1  JV“1 

N 

-  AT-1  N 

E 

X  (0<  -  00 , 

Bv  =  —  y]  hihk  (Qi^i  ~  QkVk)  • 

(A.28) 

(A.29) 

(A.30) 


i=i  jfe  =  i+i 

The  multilayer  expression  for  AT©  follows  directly  from  (A.25)  and  the  observation  that  Nr  and  A/©  differ 
ultimately  only  by  an  interchange  of  eRR  and  e%e,  and  k,RR  and  fc©©,  which  yields 

Nq  —  N  +  Ae©©  +  Av  cRR  +  B  kee  +  BvkRR. 

The  multilayer  expression  for  Nrq  is  an  obvious  generalization  of  (8.142),  i.e., 

Nrq  =  Aq  e°RQ  +  Be  Hrq  , 

where  the  constants  generalize  those  found  in  (8.98)  and  (8.99): 

A©  =  X  hiGi ,  Be  =  \Y,  X  hihk  "  G(t)  •  (A'31) 

«=i  <=  i  k  =  i+i 

Turning  now  to  the  calculation  of  Mr  for  a  multilayer,  we  have  from  (A.8): 

N  ,€i  _  *  rU  ft  "  /* 

Mr  =  X  /  (f  -  V2)  S(o)RRi  d£  =  X  /  £  5(0 )RRi  -  2  X  /  S(0)RRi  d£. 

Jt=t i_i  <=1  »=i 

Substituting  from  the  first  of  equations  (A.15),  and  using  the  definition  (A.5)  in  the  last  term,  the  last 
equation  yields 


Mr 


f  Wrn  ~  (£  -  h/2)  t]Ri ]  d£  -  -A(r, 


-tf 

=  E  f  (»*  +  s”") « -  E 

=  5E  (™  +  5>w)  («?  -  <t.)  -  jEo*  (f<  -  «■-,)  - 

or,  factoring  the  squared  and  cubed  difference  terms,  and  using  (A.2)  in  the  results: 

1  N  h  N 

Mr  =  -X  aR*  hi  (6  +  6-l)  +  4  X  ^  hi  &  +  &-0 
1  »=1  «=1 

-  $  +  m*-1  +  6-0  - 1^-  (A<32) 
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Using  the  last  equation  of  (A.-17)  to  replace  Nr  in  (A. 32),  and  collecting  terms  involving  am  and  rjm  in  that 
result,  we  obtain 


1  N 

Mr  —  -- ^2  aRi  hi  {h  ~  Ci  -  £i-i) 


i= 1 


1  N 

~  Yo  E  ^Rihi  [3h*  -  6ft  (&  +f?-i)  4*  4  (£f  4-  +  £?_i)] 

1  i=  1 


or,  using  (A.20)  and  (A.22)  in  the  sum  involving  or^ 


1  N—l  N 

Mr  =  ~2  E  E  hihk (*»  -  ***) 

»=  1  fc  =  i+l 

1  * 

“  TK  E  ^Ri  ^  i^h2  —  6ft  (£*  4-  Ci-i)  +  4  (ff  +  &6-i  4-  Cf-i)]  *  (A33) 

^  i=l 

Substituting  for  aRi  and  7}m  from  (A.13),  we  find  that  (A. 33)  can  be  written  in  the  same  form  as  equation 
(8.144),  Le., 

Mr  =  -M-  Be°RR  -  Bv €©@  -  DkRR  -  A,fe©e,  (A.34) 

where  the  multilayer  constants  B  and  Bv  are  defined  in  (A.28),  M  is  given  by 

1  Ar— 1  N 

M=2T,  E  bhk(Si-Sk),  (A.35) 

i~l  k  =  i-f  1 

and  D.  Dv  have  the  forms 

1  N 

D  =  12  E  Qi  hi  l3h2  -  6h  (6  +  6- 1)  +  4  (£?  +  66-1  +  f?_i)]  ,  (A.36) 

*=1 

i  * 

A/  =  12  E  Q'vihi  [3 ft2  -  6ft  (6  +  6-1)  +  4  (€?  +  66-1  +  Cj-i)]  ■  (A.37) 

These  two  coefficients  can  be  simplified  somewhat  using  the  following  identities: 

3 ft2  -  6ft  (6  +  6-i)  =  3  [ft2  -  2 ft  +  26-0]  s  3  [(ft  -  ft*)2  -  ft?  -  4ft6-i] , 

£<  +  66-1  +  6-1  =  (6  -  6-i)2  +  366-1  =  ft?  +  36  6-i, 

where  we  made  use  of  the  fact  that  &  +  &_i  =  hi  4-  in  the  first  of  these.  Substituting  these  results  in, 
for  example,  (A.36),  yields  after  simplifying: 


D  = 


^  E  Qi  fti  [ft?  +  3  (ft  -  hi)2  -  126-1  (ft  -  6)] , 


(A.38) 


and  we  note  that  the  third  term  in  brackets  vanishes  for  i  =  1  (since  £0  =  0)  and  i  =  N  (since  =  ft). 
Similarly,  Dv  may  be  written  as 


1  ^ 

Dv  =  i2  E  **  ^  +  3  (ft  -  fti)2  -  126-1  (ft  -  6)] 


i=l 


(A.39) 
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The  multilayer  form  for  Me  follows  directly  from  (A. 34)  by  simply  interchanging  e°RR  and  e%e,  and  kRR 
and  kee ,  which  yields 

M©  =  ~M  -  Be°ee  -  Bue°RR  -  Dkee  -  DvkRR ,  (A.40) 

while  the  multilayer  expression  for  MRq  is  a  generalization  of  (8.145),  i.e., 

MRe  =  -BeeRe  ~~  D&kRe,  (A.41) 

where  Be  is  given  in  (A. 31),  and  the  new  constant  De  generalizes  (8.112)  (it  is  in  fact  the  same  as  (A.38) 
with  Qi  replaced  by  <?*): 

1  N 

=  [h2i  +3  (h-  hi)2  -  12^_i  ( h  -  &)] .  (A.42) 

t=l 

It  is  straightforward  to  check  that  our  earlier  results  for  the  stress  resultants  and  couples  (for  a  membrane 
with  a  single  coating)  follow  immediately  from  those  determined  in  this  Appendix  by  taking  N  =  2  (in  which 
case  i  =  1  ==  c  corresponds  to  the  coating,  while  i  =  2  =  s  corresponds  to  the  membrane  substrate). 

Finally,  for  a  multilayer  coating  the  areal  density  70  defined  by  equation  (8.103),  which  determines  the 
gravitational  body  force  appearing  in  the  axial  equilibrium  equation  of  each  theory,  must  be  generalized  to 

TV 

Jo  =  ^hipoi,  (A.43) 

t=l 

where  poi  is  the  mass  density  of  the  material  in  layer  i. 
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